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CHAPTER  1 


INTRODUCTION 

The  interaction  of  gases  with  solid  surfaces  is  a  topic  of  considerable 
fundamental  as  well  as  technical  interest.  Such  processes  as  catalytic 
reactions,  corrosion,  and  oxidation  begin  with  the  collision  of  a  gas 
molecule  with  a  surface.  One  of  the  ways  to  study  this  Interaction  is  mo¬ 
lecular  beam  scattering.1'  In  this  technique,  a  well  defined  beam  of  atoms 
or  molecules  is  allowed  to  collide  with  a  crystal  surface  in  an  ultrahigh 
vacuum.  The  resulting  scattering  distribution  contains  considerable  infor¬ 
mation,  not  only  on  the  scattering  process  itself,  but  also  on  the  potential 
and  structure  at  the  surface. 

In  this  thesis,  we  investigate  several  aspects  of  the  scattering  of 
light  molecules  (H^ ,  D, ,  HD,  and  He)  at  thermal  energies  from  magnesium 
oxide.  The  magnesium  oxide  surface  is  of  interest  not  only  as  an  elementary 
oxide,  but  also  as  a  simple  model  system  on  which  to  study  the  dynamics  of 
the  gas-surface  interaction.  Our  primary  emphasis  is  on  the  elastic  and 
inelastic  diffraction  of  the  molecular  beam  from  the  crystal  surface.  Of 
particular  interest  are  those  inelastic  processes  involving  rotational 
energy,  which  can  accompany  the  scattering  of  diatomic  molecules.  We  will 
show  that  a  new  effect,  extremal  scattering,  can  cause  singularities  in  the 
intensity  of  these  inelastic  rotational  processes.  Our  chief  concern  is 
in  studying  those  scattering  processes  which  are  unique  to  molecular  species, 
that  is  processes  involving  Che  Internal  degrees  of  freedom.  We  have, 
however,  undertaken  limited  studies  of  helium  diffraction  from  magnesium 
oxide  to  afford  a  baseline  for  comparison  with  the  more  complex  molecular 
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Molecular  beam  scattering  in  the  gas  phase,  either  single  beam  or 
crossed  beam,  has  been  used  in  chemical  investigations  for  many  years. 

It  has  yielded  considerable  information  not  only  on  the  intermolecular 
potentials,  but  also  on  energy  transfer  and  the  reactivity  of  vibrational ly , 
rotationally ,  and  electronically  excited  species.  In  contrast,  scattering 
from  surfaces  has  not  yet  reached  this  level  of  sophistication.  The 
scattering  process  can  be  completely  characterized  by  the  translational 
energy,  the  internal  state,  and  the  angular  distribution  of  both  the 
incident  and  scattered  molecules.  For  gas  phase  and  beam-beam  scattering, 
these  can  be  measured  or  controlled  with  relative  ease.  For  instance, 
lasers  can  be  used  to  selectively  excite  the  internal  states  of  the  beam, 
and  fluorescence  used  to  monitor  the  state  of  the.  scattered  species. 
Unfortunately,  these  techniques  cannot  in  general  be  applied  in  studying 
scattering  from  surfaces.  The  requirements  of  ultrahigh  vacuum  keep  the 
beam  and  scattered  intensities  far  too  low.  At  present,  for  the  most  part, 
we  must  be  content  with  knowledge  of  the  average  incident  energy  and  angle, 
and  the  scattered  angles.  We  must  to  a  large  extent  rely  on  interpretation 
of  peak  shapes  and  positions  in  the  scattered  distribution  to  arrive  at  in* 
sights  about  surface  Interactions.  This  restriction  severely  limits  the 
type  of  information  available  on  inelastic  processes  at  surfaces;  for  the 
most  part,  only  qualitative  information  is  available  on  these  processes. 

Diffraction  of  a  molecular  beam  was  first  observed  by  Stern  and 

Estermann,^  some  40  years  ago,  for  H,  and  helium  from  lithium  fluoride. 

This  work  confirmed  the  deBroglie  hypothesis  for  atoms,  coming  only  shortly 
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after  the  diffraction  of  electrons  was  observed  by  Davisson  and  Germer. 


The  experiments  of  Stern  and  co-workers  also  yielded  information  on  gas- 
surface  interactions,  and  provided  a  considerable  stimulus  for  the  early 
theories  of  surface  behavior. 

Despite  the  success  of  Stem  and  Estermann,  relatively  little  work 
on  the  dynamics  of  gas-surface  interactions  was  done  using  scattering 
methods.  The  development  of  molecular  beam  scattering  as  a  surface 
technique  had  to  await  advances  in  vacuum  and  electronic  technology.  The 
60' s  marked  the  beginning  of  modern  work  on  the  scattering  of  rare  gases 
as  well  as  hydrogenlc  molecules  from  both  ionic**’^’6  and  metal^’6  surfaces. 

Much  of  this  early  work  was  concentrated  on  the  diffraction  of 
helium-4’6  from  the  alkali  halide  surfaces.  However,  diffraction  of  hydro- 
genic  molecules'*’6  was  also  observed.  In  general,  the  angular  distributions 
for  D,  were  found  to  be  broader  than  those  for  This  was  believed  to 

be  the  result  of  inelastic  processes  between  rotational  modes  of  the 
molecules  and  surface  phonons.  The  rotational  spacing  of  is  approxi¬ 
mately  one-half  that  of  end  is  nearer  the  Debye  phonon  energy  of  these 
surfaces.  Because  of  the  continuous  energy  range  of  the  phonons,  it  was 
proposed  the  such  inelastic  processes  would  result  in  a  smearing  of  the 
scattering  features.  There  was,  however,  no  direct  evidence  for  this  type 
of  rotational  energy  exchange. 

9 

Quite  another  type  of  rotational  energy  transfer  is  possible.  Logan 
pointed  out  some  time  ago  that  rotational  and  translational  energy  could 
be  directly  Interchanged  in  scattering  at  a  surface.  He  noted  that  if  such 
a  process  occurred,  it  would  result  in  discrete  peaks,  shifted  from  the 
specular,  containing  only  molecules  of  a  given  rotational  state.  The 
possibility  of  inelastic  diffraction  was  not  considered,  however;  only 
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transitions  from  the  specular  peak  were  considered.  Little  notice  was  taken 
of  this  work,  as  these  peaks  were  not  observed  in  the  data  available  at  the 
t  iae . 

Direct  evidence  for  this  type  of  process  in  gas-surface  scattering  was 

first  obtained  in  1974  by  R.  Grant  Rowe*0, at  Illinois  and  independently 
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in  Genoa  by  Soata,  et  al.  Rowe  observed  it  for  H^.  D., ,  HD  from  MgO  (001). 
In  Genoa,  H-,  from  cold  (80  K)  LiF  was  studied.  Distinct  additional  peaks 
in  the  scattering  distribution  were  seen.  Both  experiments  can  be  fit 
well  to  a  model  which  postulates  energy  transfer  between  rotational  and 
trans lational  modes,  without  transfer  to  the  surface.  In  simple  terms, 
this  model  can  be  viewed  as  scattering  of  a  rigid  rotor  from  a  rigid  periodic 
surface.  The  phenomenon  is  known  as  rotational  inelastic  diffraction.  The 
magnitude  of  the  rotational  transfer,  as  well  as  the  sharpness  of  the 
resulting  peaks,  were  striking  in  the  early  data  on  MgO  (001)  for  incidence 
along  the  (100]. The  experiments  in  this  thesis  were  undertaken  to 
extend  these  measurements  to  the  (110]  azimuth,  and  in  particular  to  in¬ 
vestigate  the  mechenism  responsible  for  the  sharpness  of  the  rotational 
transitions.  Singularities  in  the  scattering  probability  have  been  found 
and  explained.  As  a  result,  the  peak  shapes  and  positions  for  rotatlonally 
inelastic  diffraction  can  now  be  accurately  predicted.  All  our  experimental 
work  has  been  done  on  MgO  (001).  This  surface  is  structurally  similar  to  the 
LIF  (001),  which  is  more  commonly  used  in  bean  scattering  studies.  We  have 
found,  however,  that  due  to  its  higher  Debye  temperature,  thermal  inelastic 
scattering  is  considerably  reduced.  The  diffraction  features  are  thus  more 
clearly  defined. 

Our  primary  emphasis  is  on  the  participation  of  rotational  degrees  of 
freedom  in  the  scattering  of  molecules  from  a  surface.  We  have  also  under¬ 
taken  studies  of  helium  scattering  from  these  same  surfaces,  for  comparison 


with  the  more  complex  molecular  scattering  data.  These  data  for  helium, 
however,  are  not  as  simple  to  Interpret  as  originally  expected.  We  have 
therefore  carried  out  calculations  of  the  diffraction  profiles  using  the 
seaic lasslcai  approximation. 

Although  several  quancum  mechanical  formalisms^  ^  for  gas-surface 
scattering  have  been  proposed,  most  are  cumbersome,  If  not  Impossible  to 
use  In  realistic  situations.  The  numerical  calculations  are  extensive, 
and  convey  little  feeling  for  the  mechanics  of  the  Interaction.  Theoretical 
chemists  have  In  the  past  used  semlc lass  leal  calculations  for  beam-beam 
scattering  at  low  energies.^'  The  semlc lass  leal  formalism  Is  an  outgrowth 
of  Che  classic  UKB  method  and  has  been  extensively  developed  for  scattering 
scudles  by  Marcus*  and  by  Miller.  In  many  cases.  It  gives  quantitatively 
accurate  results  for  Intensities  In  gas-surface  scattering. The  calcu¬ 
lations  are  not  difficult  and,  as  a  result,  the  formalism  Is  Ideal  for 

fitting  and  Interpreting  experimental  data.  Our  calculations  will  be  done 
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using  a  simple  model  potential,  known  as  Steele's  potential.  This  model 
potential  Is  sufficient,  however,  to  allow  identification  of  the  essential 
teatures  in  Che  helium  scattering  data. 

Although  in  the  past,  work  in  Che  area  of  gas-surface  scattering  was 
limited,  there  has  been  considerable  activity  in  the  last  5  years.  A  new 
generation  of  very  sophisticated  molecular  beam  scattering  systems  has  been 
developed,  with  instrumentation  and  capabilities  far  surpassing  those  of  our 
present  apparatus.  Thus,  after  these  initial  studies  on  MgO,  we  have 
modified  our  system  for  work  in  s  new,  relatively  unexplored  area.  We  have 
developed  a  method  of  using  the  unique  capabilities  of  the  molecular  beae. 
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to  study  surface  diffusion  of  gaseous  species  on  crystal  surfaces,  about 
which  only  little  is  known.  This  new  technique  will  be  discussed  briefly 
in  Chapter  6. 
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CHAPTER  2 
FUNDAMENTALS 

2.1  Surface  Geometry 

Magnesium  oxide  is  an  ionic  crystal  having  the  NaCl  structure.  It 
cleaves  easily  along  the  (001)  plane,  the  structure  of  which  is  shown  in 
Fig.  1.  The  (001)  is  a  non-polar  plane,  having  equal  numbers  of  positive 
and  negative  ions,  because  of  the  discrete  ion  cores,  the  potential  varies 
strongly  across  the  surface.  This  periodic  variation  in  the  potential  is 
responsible  for  the  relatively  strong  diffraction  observed  from  these 
surfaces.  The  surface  unit  cell,  shown  in  Fig.  1,  is  2.97  A  on  a  side 
and  is  rotated  43*  from  the  bulk  unit  cell.  The  surface  reciprocal  lattice 
vectors,  also  shown  in  Fig.  1,  thus  lie  along  the  (110]  directions.  In 
these  studies,  we  have,  for  the  most  part,  investigated  scattering 
along  the  [110]  incident  azimuth.  In  other  words,  the  surface  component 
of  the  incident  wavevector  lies  along  the  [110]  direction,  the  direction 
between  nearest  like  neighbors. 

In  scattering  studies,  the  Incident  plane  is  that  defined  by  the  incident 
wave  vector  and  the  surface  normal;  the  scattered  plane  is  that  defined  by 
the  incident  wave  vector  and  the  scattered  wave  vector.  The  simplest  case 
is  that  in  which  the  incident  and  scattered  planes  coincide;  this  is  known 
as  in-plane  scattering.  Although  the  beam  may  be  scattered  into  any  portion 
of  the  half-space  above  the  crystal,  considerable  information  can  be  gained 
by  considering  only  the  ln-plane  portion.  The  geometry  appropriate  to  these 
studies  of  ln-plane  scattering  is  shown  in  Fig.  2,  for  the  particular  case 
of  incidence  along  the  [110]  azimuth.  The  incident  and  scattered  angles, 
measured  up  from  the  surface,  are  denoted  by  9Q  and  9  respectively;  the 
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experimentally  measured  angle  T*t,  the  angle  between  the  incident  and 
scattered  beams,  is  therefore  given  by 

-  180  -  90  -  e.  (1) 


2.2  Diffraction 

In  our  scattering  experiments,  we  are  for  the  most  part  concerned  with 

the  diffraction  of  the  molecular  beam  from  the  periodic  surface  potential. 

According  to  the  deBroglie  hypothesis,  the  particles  will  have  a  wavelength 

X.  -  h/p.  For  He  atoms  near  room  temperature  the  most  probable  wavelength, 

X  ,  is  thus  approximately  0.5  A.  Since  the  surface  lattice  has  a  periodicity 

np 

on  the  order  of  3  A,  we  can  expect  diffraction  of  the  molecular  beam  from 
the  surface,  if  the  scattering  is  sufficiently  coherent.  It  should  be  noted 
that  at  thermal  energies,  that  is,  translational  energies  in  the  range  of 
40  meV,  scattering  occurs  from  the  repulsive  hard  wall  potential  of  the 
surface  layer  of  atoms;  the  incident  atoms  do  not  penetrate  the  surface. 


2.3  Diffraction  of  Monoenerget lc  Beam 

The  scattering  of  a  monoenergetic  beam  is  simplest  and  will  be  dis¬ 
cussed  first.  We  may  describe  the  incident  and  scattered  beams  by  their 

*4  -* 

wave  vectors,  k^  and  k,  respectively.  The  surface  components  of  these  are 

denoted  by  K  and  3,  that  is 
o 


-  (l.s  z)  . 

o  o  oz 


The  Laue  condition  for  a  2-D  potential  Is  thus 


(2) 


4  3 


nm 


(3) 


where  is  a  surface  reciprocal  lattice  vector.  For  a  square  symmetric 
-♦  2tt  «  2n 

lattice,  G _  ■  -T-  nx  4  —  my.  This  condition,  along  with  conservation 

ran  a  a 
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2  2 

of  energy  k^  •  k  ,  uniquely  determines  the  positions  of  the  diffracted 
peaks.  For  the  particular  case  of  in-plane  scattering,  we  have 

cos  9  ■  cos  9  ♦  -r  •  (4) 

o  d 

The  scattering  distribution  thus  consists  of  a  specular  peak,  that  is  the 
(0,0)  diffraction  peak,  surrounded  by  a  series  of  discrete  diffraction 
peaks.  This  is  the  situation  appropriate  to  the  scattering  of  low  energy 
helium  atoms  from  an  ionic  surface. 


2.4  Inelastic  Diffraction 

In  addition  to  elastic  scattering,  several  types  of  inelastic  events 

may  occur.  A  molecule  can  exchange  energy  from  either  translational  or 

internal  modes  with  the  surface  by  coupling  to  phonons.  At  the  energies 

we  are  considering,  electronic  and  vibrational  excitations  are  forbidden; 

however,  rotational  degrees  of  freedom,  having  energies  in  the  range  10 

to  100  meV,  can  freely  participate.  Another  type  of  inelastic  event  Involves 

direct  transfer  of  energy  between  translation  and  internal  modes,  without 

energy  exchenge  with  the  surface.  If  this  occurs  during  diffraction,  we 

10-12 

have  rotationally  inelastic  diffraction.  For  this  we  have 


k 


2 


2M  A E 


(5a) 


(5b) 


where  AE  is  the  energy  transferred  from  rotation  to  translation.  Since  the 
rotational  energy  levels  are  quantized,  new  discrete  Inelastic  peaks  are 
expected.  Transitions  in  which  energy  is  transferred  from  rotation  to 
translation  (AE  >  0)  are  called  gain  transitions;  if  energy  is  transferred 
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from  translation  (&£  <  0)  they  are  called  loss  transitions.  Confining 
ourselves  again  to  in-plane  scattering,  the  peak  positions  are  given  by 

cos  9  ■  (cos  8  +  nX  /d)  (6a) 

^  oo 

o 

where 

X A  -  (1  2&Z  M  X2/h2]‘5  .  (6b) 

o  *  go 

For  the  symmetric  rotors  H,  and  ,  symmetry  forbids  transitions 

between  even  and  odd  rotational  states,  i.e.,  there  is  no  ortho-para 
25 

conversion.  The  lowest  allowed  transitions  are  those  in  which  the 
rotational  quantum  number  J  changes  from  0  -♦  2 ,  followed  by  l-»3,  then 
2  -*  l*.  For  HD,  however,  there  is  no  such  selction  rule.  Transitions 
between  adjacent  rotational  states  are  allowed.  The  spacings  of  Che 
rotational  transitions  for  the  hydrogenic  molecules  are  given  in  Tables 
1-3,  along  with  the  equilibrium  population  P(J^),  of  the  initial  state, 

J^,  at  two  beam  temperatures,  300  K  and  80  K.  The  transitions  are  labeled 
by  1*4-1,  ±2,  ♦  3  for  gain  (+)  and  loss(-)  transitions,  the  integer  denoting 
Che  relative  magnitude  of  the  energy  change  in  the  transition.  An  inelastic 
diffraction  peak  is  thus  fully  described  by  (n,m)  +  I.  The  energy  change 
in  the  lowest  allowed  transition  for  HD  is  seen  to  be  255  cal/mode,  compared 
to  512  cal/mole  and  1012  cal/tnole  for  D,  and  H.,  respectively.  By  comparison, 
Che  mean  translational  energy  of  these  molecules  in  a  beam  at  300  K  is 
1190  cal/mole.  Considerable  changes  in  Che  positions  of  the  diffracted 
peaks  can  therefore  be  expected  due  to  these  rotational  transitions. 

2.5  Diffraction  of  a  Thermal  Beam 

We  have  discussed  the  diffraction  of  a  monoenergetlc  beam  of  molecules. 
All  beams,  however,  have  a  distribution  of  velocities.  A  fairly  mono- 
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TABLE  1 

Rotational  energies  and  allowed  rotational  transitions  for  HD. 


Rotational  Energy  Levels  (HD) 

J 

0 

1 

2 

3 

4 

Allowed  Transitions  (HD) 


ROT 

0.0  cal/mole 
255.06 
763  .Ui* 

1521.65 

2524.6 


I 

VJf 

4E(I) 

cal/mole 

PU*) 

300  K 

80  K 

+1 

1  -  0 

+  255.06 

0.39 

0.37 

+2 

2  -  l 

a-  508.38 

0.28 

0.025 

+3 

3  -  2 

+  758.21 

0.11 

— 

-l 

0-  1 

-  255.06 

0.20 

0.61 

-2 

1-2 

-  508.38 

0.39 

0.37 

-3 

2-3 

-  758.21 

0.27 

0.025 

-4 

3-4 

-1002.95 

0.11 

0.025 
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TABLE  2 

Rotational  energies  and  allowed  rotational  transitions  for  D2* 

Rotational  Energy  Levels  (H^) 

J  E_  „ 

ROT 

0  0.0  cal/mole 

1  338.68 

2  1012.93 

3  2016.64 

4  3340.39 


Allowed  Transitions  (H,) 


I 

Ji*Jf 

AE 

cal/mole 

P(JL) 

300  K 

80  K 

+1 

2-0 

+  1012 

0.117 

0.002 

♦2 

3-  l 

+  16  78 

0.091 

— 

+3 

4-2 

+  2327 

0.004 

— 

-1 

N 

t 

o 

-  1012 

0.128 

0.248 

-2 

1-3 

-  16  78 

0.658 

0.750 

-3 

2-4 

-  2327 

0.117 

0.002 

TABLE  3 

Rotational  energies  and  allowed  rotational  transitions  for  H. 


Rotational  Energy  Levels 

J 

0 

1 

2 

3 

4 


0>2> 

erot 

0.0  cal/mole 
170.8 
511.7 
102  1 . 1 
1696.6 


Allowed  Transitions  (D,) 

I  jrjt 

AE(I) 

cal/mole 

300  K 

+1 

2  -*0 

♦512 

0.38 

♦2 

3  —  l 

♦850 

0.11 

+3 

4  -*  2 

♦  1185 

0.09 

-l 

0-2 

-512 

0.181 

-2 

1-3 

-850 

0.204 

-3 

2-4 

-  1185 

0.38 

P(JL) 


80  K 
0.111 
0.003 

0.554 
0.329 
0.  Ill 


l 


0.003 
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energetic  beam  of  neutrals  can  be  obtained  by  aerodynamic  expansion  from 

26  2  7 

high  pressure,  a  so-called  nozzle  beam.  '  With  such  a  source,  velocity 

half-widths  on  the  order  of  1%  can  be  obtained  from  room  temperature  beams. 

A  much  broader  velocity  distribution  is  obtained  in  a  beam  formed  by 
28  29 

effusion."  ’  Molecules  are  allowed  to  effuse  from  a  low  pressure  source 
into  vacuum  through  a  small  hole  or  channel.  If  the  mean  free  path  in  the 
gas  behind  the  source  is  much  longer  than  the  diameter  of  the  hole,  the 
distribution  of  velocities  in  the  beam  will  reflect  the  equilibrium 
velocity  distribution  of  gas  in  the  source.  The  beam  will  have  a  modified 
Maxwe 11-Boltzman  distribution,  given  by 

Jt  2  3  V! 

g(v)dv  -  2(^£*)  v  exp (-  .  (7) 

where  T  is  the  temperature  of  the  source.  The  mean  translational  energy 
of  the  beam  molecules  is  2kT.  Recalling  that  X.  •  h/M^v,  the  wavelengths 
in  the  beam  are  distributed  according  to 


f  (X)dX  -  -2  (X  /\)5exp(-V*/\2)  d  (~)  , 

Ot 


(8) 


2  h 

where  *  77—  .  The  most  probable  wavelength  \  is  thus  given  bv 

Or  L  KTM  mp 

s 


(9) 


The  diffraction  of  a  molecular  beam  from  the  2-D  surface  potential  is 
in  many  ways  analogous  Co  Che  diffraction  of  X-rays  from  the  bulk  crystal. 
For  the  latter  we  have  the  3-D  Laue  condition  for  constructive  inter¬ 


ference  of  the  scattered  waves. 
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k  -  ic  +  2  .  (io) 

o  run 

For  a  given  sec  of  incidenc  angles,  only  certain  wavelengths  will  diffract; 
all  others  destructively  interfere.  Because  of  this,  one  can  do  X-ray 
diffraction  from  a  crystal  with  a  broad  wavelength  source  and  still  get 
sharp  diffraction  spots.  The  crystal  filters  out  the  other  wavelengths. 

The  diffraction  condition  in  2-D  is  much  less  strict.  Only  the  two 

•4  *4 

equations  K  •  and  conservation  of  energy  must  be  satisfied.  As 

a  result,  all  wavelengths  will  be  diffracted  from  the  surface,  but  each 
will  be  scattered  to  a  different  angle.  For  a  beam  with  a  narrow  velocity 
distribution,  this  presents  no  problem.  For  a  beam  with  a  significant 
wavelength  distribution,  as  in  £q.  (8),  the  situation  is  more  complex. 

Broad  diffraction  peaks  are  then  expected,  reflecting  the  wavelength 
distribution  of  the  incident  beam.  This  is  shown  schematically  in  Fig.  3, 
where  calculated  profiles  for  typical  (+  1,0)  diffraction  peaks  are  shown, 
assuming  that  the  diffraction  probability  is  independent  of  wavelength. 

It  should  be  noted  that  even  in  this  simple  case,  the  peak  does  not 
correspond  to  the  most  probable  wavelength,  or  to  any  other  simple 
characteristic  of  the  beam.  The  wavelength  distribution,  the  transition 
probability,  and  a  geometrical  factor,  d\/d@,  all  interact  to  determine  the 
peak  of  the  scattered  distribution.  The  situation  is  even  more  complex  for 
inelastic  scattering.  It  it  precisely  this  geometrical  factor  which  will 
have  a  significant  effect  on  the  detectability  of  the  rotational  inelastic 
events  we  wish  to  study.  This  will  be  fully  investigated  in  Chapter  5. 


Simulated 


Scattering  Angle,  17  (degrees) 


Calculated  first  order  helium  diffraction  peaks  for  thermal  velocity 
distribution  at  300  K,  assuming  constant  diffraction  probability. 

The  broad  peaks  reflect  the  broad  initial  wavelength  distribution. 
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CHAPTER  3 

EXPERIMENTAL  APPARATUS  AND  TECHNIQUE 

Molecular  bean  experiments  are  conceptually  simple.  ^  One 
merely  scatters  a  beam  from  a  surface  and  measures  the  scattered  Intensity 
as  a  function  of  angle.  In  practice,  however,  considerable  difficulties 
arise.  To  do  molecular  beam  scattering,  and  in  particular  to  study 
diffraction,  one  needs  a  well  defined,  preferably  coherent,  source  of 
low  energy  molecules.  An  ultrahigh  vacuum  chamber  is  necessary,  In  which 
a  clean  sample  can  be  prepared  and  maintained.  One  also  needs  an 
efficient  method  of  detecting  the  scattered  neutral  molecules  and  discrim¬ 
inating  against  the  background.  The  apparatus  described  below  accomplishes 
these  objectives. 

Our  molecular  beam  system  consists  of  3  sections;  the  main 

scattering  chamber,  the  beam  formation  chamber,  and  the  gas  handling  section. 

The  entire  apparatus  is  constructed  of  stainless  steel,  OFHC  copper, 

o 

and  pyrex  glass,  and  can  be  routinely  baked  to  300  C.  Ultrahigh  vacuum  is 
attained  In  all  sections.  Schematic  views  of  the  system  are  shown  in 
Figs.  4  and  5. 

3. 1  Gas  Handling  Svstem 

To  operate  a  molecular  beam,  gas  must  be  delivered  to  the  source 
at  a  constant  pressure.  Since  the  beam  flux  hitting  the  sample  is  typically 
a  mono  layer/ sec ,  even  small  amounts  of  impurities  could  cause  contamination. 

A  system  Is  therefore  required  both  to  purify  the  gas  and  establish  a 
constant  pressure.  A  schematic  of  the  apparatus  built  for  this  purpose  Is 
shown  in  Fig.  6.  This  section  is  constructed  entirely  of  pyrex,  with 
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Fig.  6.  Diagram  of  gas  handling  system  with  automatic  pressure  control  (APC) 

p  and  P  ,  are  the  measuring  and  reference  (dummy)  Pirani  gauges 
meas  ref 

respectively.  Gas  is  stored  in  the  reservoir,  R. 


metal  sealed  L'HV  values  (Granville-Philllps,  type  C  and  series  204)  and  Is 
evacuated  by  a  mercury  diffusion  pump  ( Konces-Martin ) .  DeuCerlum  and 
helium  are  drawn  from  analyzed  reagent  grade  flasks  (Alrco).  H,,  is  obtained 
by  diffusion  through  a  palladium  leak  (K  6  B  Apparatus,  Schenectady,  NY). 

HD  is  made  by  equilibracing  a  50/50  Hj-D^  mixCure  over  *  hot  tungsten 
filament,  resulting  in  a  mixture  of  H^ ,  D^,  and  HD.  These  gases  are  stored 
in  the  reservoir  at  a  pressure  of  several  hundred  torr. 

For  proper  operation,  the  source  requires  a  constant  pressure  of 
about  1  torr.  It  is  necessary  to  control  this  pressure  over  periods  of 
hours  to  +  l!i.  This  is  accomplished  by  a  servo  valve  (Cranville-Phillips , 
Series  216)  using  two  Pirani  gauges  (LKB,  3596)  as  a  differential  pressure 
sensor.  One  Pirani  gauge  is  blanked  off  at  the  desired  pressure.  The 
other  measures  the  pressure  delivered  to  the  source.  The  outputs  of  the 
two  are  compared  using  the  bridge  circuit  shown  in  Fig.  7.  This  arrange¬ 
ment  is  necessary  to  compensate  for  fluctuations  in  room  temperature, 
which  severely  affect  the  reading  of  any  heat  conductivity  gauge,  such  as 
a  Pirani. 


3.2  3eam  Source  and  Beam  Chamber 

The  moleculjr  beam  19  formed  by  effusion  from  a  multicapillary 

array  (Brunswick  Cor p. ,  Collimated  Hole  Structure).  To  form  a  beam  by 

effusion,  it  is  necessary  to  operate  in  the  region  where  the  mean  free 

path,  L  .  ,  is  greater  than  d,  the  diameter  of  the  channel.  Because  of 
mtp 

this  restriction,  more  intensity  can  be  achieved  by  using  an  array  of 
small  channels  rather  than  a  large  channel.  This  is  the  rationale  behind 
the  multicapillary  source. ^  The  source  consists  of  about  100  parallel 


Ov’pv'  *9* 

Co«"*o»  C*»c-  * 


Schematic  diagram  of  differential  Pirani  bridge  circuit 
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channels,  each  3  am  long  and  . 05  mm  In  diameter,  through  a  stainless  steel 
disc.  The  channels  are  close-packed  In  a  circular  array  about  1  mm  In 
diameter.  A  diagram  of  the  source  assembly  Is  shown  In  Fig.  8.  The 
source  disc  Is  brazed  In  the  end  of  the  source  tube,  forming  a  small 
effusion  chamber.  Gas  Is  supplied  to  this  chamber  at  constant  pressure 
by  the  gas  handling  system  previously  described.  The  source  tube  can  be 
cooled  to  change  the  velocity  distribution  of  the  beam.  This  Is 
accomplished  by  flowing  cold  nitrogen  gas  through  the  colls  around  the 
source  tube.  A  window  Is  provided  at  the  rear  of  the  source  tube  for 
alignment  of  the  source  and  carget  with  a  laser. 

The  molecular  beam  for  our  experiments  Is  formed  by  gas  effusing 
through  Che  source  channels  Into  the  vacuum.  It  thus  has  a  velocity 
distribution  given  by  Eqn.  (7)  and  Is  peaked  In  the  forward  direction. 

A  collimator  hole,  0.7  ms  In  diameter,  approximately  1  cm  downstream  from 
Che  source,  allows  only  a  small  portion  of  these  molecules  to  enter  the 
main  scattering  chamber.  The  rest  of  Che  gas  Is  rejected  back  Into  the 
beam  chamber,  where  It  Is  removed  by  differential  pumping  with  an  ultra- 
high  vacuum  diffusion  pump  stack,  (Edwards  High  Vacuum,  ITHVM2A  and  £M2). 
When  the  beam  Is  on,  Che  pressure  In  the  beam  chamber  Is  typically 
10  6  torr.  Typically,  the  beam  flux  at  the  sample  Is  10^/sec  into  a 
spot  1  mm  In  diameter. 

As  previously  discussed,  this  beam  has  a  broad  velocity 

26  2  7 

distribution.  The  alternative  source,  an  aerodynamic  nozzle,  ’  is 
characterized  by  a  much  narrower  velocity  distribution.  However,  this 
source  Is  somewhat  more  difficult  to  fabricate  and  operate  chan  a  conven¬ 
tional  effusion  source;  moreover,  the  beam  Is  more  prone  to  contamination. 


Fig.  8.  Schematic  diagram  of  molecular  beam  source  assembly 
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due  to  che  high  pressure  (~  2-10  acm)  ac  which  ic  is  operaced.  Several 

scages  of  high  speed  differencial  puraping  would  be  needed  Co  handle  che 

large  gas  loads  which  accompany  chis  source.  Furthermore,  since  we  are 

inceresced  in  scudying  che  parcicipacion  of  che  incernal  degrees  of 

freedom  of  molecules  in  che  scaccering  process,  Che  nozzle  sources  have 

a  serious  disadvancage .  The  rapid  expansion  inco  vacuum  from  high 

pressure  lowers  che  random  velocities,  chac  is  che  crans lacional 

cemperacure,  of  che  molecules.  This  energy  goes  inco  uniform  cranslacion 

of  che  bulk  gas,  resulcing  in  a  fairly  monoenergedc  beam.  In  che  same 

manner,  che  incernal  degrees  of  freedom  are  co  some  excenc  frozen  ouc  in 

32  33 

che  expansion.  The  resulcing  incernal  energy  discribucion  is  cold  ’ 
and  is  noc  well  characcerized .  In  concrasc,  che  incernal  energy  discribu¬ 
cion  of  an  effusive  beam  is  well  escablished  at  all  cemperacures .  For 
chese  reasons,  chere  appeared  co  be  some  advantage  co  using  an  effusive 
source  for  chese  scudles. 

In  order  to  facilitate  signal  processing,  the  beam  is  chopped  by 

a  25  Hz  sagneCic  tuning  fork  (American  Time  Produces,  L40C)  as  ic  enters 

che  main  chamber.  This  allows  che  use  of  lock-in  detection  of  che 

34 

scattered  signal.  The  tuned  feedback  amplifier,  shown  in  Fig.  9, 
drives  che  chopper  at  ics  resonant  frequency  and  provides  a  reference 
signal  for  che  lock-in  amplifier. 

3.3  Main  Chamber 

The  scaCCering  chamber  consists  of  a  stainless  steel  bell  jar 
30  cm  in  dlamecer.  Ic  is  evacuated  by  two  separate  pumping  systems. 

Main  pumping  is  by  a  mercury  diffusion  pump  stack  (Edwards  High  Vacuum, 
UHVM2A  and  EM2 ,  having  a  pumping  speed  of  70  liters/sec  for  air. 
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Fig.  9.  Schematic  diagram  of  drive  circuit  for  tuning  fork  chopper. 
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When  Che  molecular  beam  is  turned  on,  additional  pumping  is  required  to 

maintain  a  good  vacuum  and  a  clean  sample.  For  this  reason,  a  liquid 

helium  cryo-sorptlon  pump  (Excallbur  CW-2010)  is  built  into  the  bottom 

of  the  scattering  chamber.  When  cooled  to  4  K,  it  provides  an  additional 

pumping  speed  of  2000  llcers/sec  for  helium.  With  this  pump  activated,  the 

-10  -9 

base  pressure  is  below  10  corr,  and  remains  below  10  torr  with  the 
beam  on. 

The  molecules  scattered  from  the  crystal  are  detected  by  a 
quadrupole  mass  spectrometer  (Extranuclear  324*9).  This  allows  scattering 
experiments  with  mixed  isotopes  such  as  HD,  which  cannot  be  prepared  in 
pure  form.  The  mass  spectrometer  rotates  in  an  arc  around  the  sample. 

The  entrance  aperture  of  this  detector  determines  the  angular  resulutlon 
of  our  experiments.  This  aperture  may  be  adjusted  from  the  outside 
during  a  run  to  be  either  1.6  or  3.25  degrees.  The  rotation  of  the 
decector  is  controlled,  from  the  outside  ,  by  a  stepping  motor  through  a 
bellows  sealed  rotary  feedthrough.  Since  the  quadrupole  mass  filter 
can  only  analyze  ions,  the  neutral  molecules  from  the  beam  must  first 
be  ionized.  This  is  accomplished  with  an  electron  impact  ionizer  atop 
the  mass  filter.  Molecules  in  the  ionization  region  are  ionized  by  a 
beam  of  low  energy  electrons  and  then  electrostatically  focused  into  the 
quadrupole  filter  section.  A  21  stage  Cu-Be  electron  multiplier 
(EMI,  Ltd.)  at  the  exit  of  the  mass  filter  amplifies  the  ion  signal  by 
a  factor  of  approximately  10^  before  it  is  fed  through  the  vacuum 
enve lope . 
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3.4  Crystal  Holder  and  Specimens 

The  sample  holder  is  shown  in  Fig.  10.  The  crystal  (A)  is 
mounted  in  a  square  hole  within  a  stainless  steel  block  (B>  in  which  it 
is  free  to  slide.  It  is  held  in  place  by  small  spring  loaded  stainless 
steel  bails.  This  allows  the  crystal  to  be  moved  forward  in  the  block 
for  multiple  cleavages.  Translation  is  accomplished  by  pushing  with 
the  cleaver  rod  (C).  The  entire  crystal  block  is  mounted  inside  a  ball 
bearing  (D),  with  the  face  of  the  crystal  flush  with  the  surface  of  the 
race.  This  allows  the  crystal  to  be  rotated  about  the  surface  normal, 
so  chat  the  incident  azimuth  can  be  changed.  This  degree  of  freedom, 
however,  cannot  be  accurately  measured.  Stops  are  therefore  provided  so 
that  the  IllO]  and  [100]  azimuths  are  reproducible. 

The  entire  assembly  is  again  mounted  on  two  small  ball  bearings 
(E),  to  allow  rotation  about  an  axis  in  Che  plane  of  the  sample.  In 
this  manner,  9  Q,  the  incident  angle,  can  be  adjusted  by  pushing  Che 
cleaver  rod.  The  assembly  is  held  at  a  fixed  incident  angle  by  pressure 
from  small  spring  loaded  balls. 

The  crystal  holder  assembly  is  mounted  on  a  bellows  sealed 
rotary  feedthrough  (F),  coaxial  with  the  beam.  By  rotation  of  the  sample 
holder  about  Che  beam  axis,  the  surface  normal  can  be  rotated  out  of  the 
horizontal  plane.  This  is  essentially  a  rotation,  about  Che  beam,  of 
the  crystal  coordinate  system  with  respect  to  the  laboratory  system. 

This  degree  of  rotational  freedom  allows  us  to  bring  any  scattered 
vector  into  the  horizontal  plane  without  changing  the  incident  conditions 
in  Che  crystal  coordinate  system.  In  this  manner,  much  of  the  scattered 
solid  angle  can  be  scanned,  while  avoiding  the  difficult  cask  of  moving 
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Che  plane  of  decaccor  rotation. 

For  meaningful  results.  It  Is  necessary  to  do  experiments  on 
a  clean  surface.  It  has  been  standard  procedure  In  beam  scattering 
experiments  to  cleave  the  crystal  external  to  the  system,4  6»^5,36 
and  lightly  heat  It  under  vacuum  to  obtain  a  clean  surface.  This  pro¬ 
cedure  Is  at  best  shaky.  The  true  state  of  such  surfaces  Is  still  open 

37  38 

to  considerable  controversy.  *  Furthermore,  thero  Is  recent  evidence 

that  the  MgO  (001)  surface  undergoes  an  Irreversible  reconstruction  upon 
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heating  to  greater  chan  4S0  K.  To  avoid  such  problems  the  crystal 

In  our  experiments  Is  cleaved  In  situ  along  the  (001)  plane.  Although 
this  Is  Che  preferred  and  rigorously  clean  approach.  It  severely  limits 
flexibility.  Heating  and  cooling  the  crystal  are  for  the  most  part 
precluded.  As  a  result,  these  experiments  are  conducted  with  the  surface 
at  room  temperature,  resulting  In  an  unfortunately  large  background  of 
thermal  Inelastic  events. 

The  magnesium  oxide  specimens  were  obtained  from  the  Norton  Co. , 
In  Che  form  of  l  x  1  x  5  cm  blocks.  The  sides  were  ground  down  to  uniform 
size,  using  silicon  carbide  paper,  to  permit  smooch  sliding  In  Che 
multiple  cleavage  sample  holder.  The  samples  were  then  thoroughly  washed 
In  ethyl  alcohol  and  dried  In  flowing  nitrogen,  prior  to  Insertion  In 
the  vacuum  system. 

3.5  Data  Handling 

The  scattered  beam  Is  detected  synchronously  with  the  chopped 
Incident  beam  using  a  lock- In  amplifier  (PAR  186).  However,  even  with  a 
time  constant  of  several  seconds  on  the  output  filter,  the  signal  is 
noisy.  Since  signal  Is  accumulated  at  each  angle  for  10-30  seconds,  it 


would  be  desirable  to  integrate  the  signal  over  this  entire  time  Interval 


to  reduce  the  noise.  To  accomplish  this,  and  to  control  the  experiment, 
a  signal  integrator  and  digitizer  was  constructed.  A  function  block 
diagram  of  chis  instrument  is  shown  in  Fig.  11.  The  main  section  of  this 
instrument  is  a  gated  integrator.  After  the  detector  is  in  its  desired 
position,  the  integrator  is  enabled.  It  integrates  the  output  of  the 
lock-in  for  a  preset  interval  while  the  detector  sits  at  one  position. 

This  period  may  be  adjusted  from  5-60  seconds.  At  the  end  of  an  interval, 
the  input  to  the  integrator  is  disabled,  and  the  accumulated  signal  is 
digitized.  Simultaneously,  the  angular  position  of  the  detector,  as 
measured  with  a  potentiometer ,  is  digitized.  Both  pieces  of  Information 
are  then  printed  onto  a  teletype  and  punched  on  paper  tape  for  later 
computer  analysis.  Finally,  the  integrator  is  reset  to  zero  and  the 
detector  is  automatically  stepped  0.45°.  The  sequence  then  begins  again 
with  the  input  of  the  integrator  being  enabled.  This  sequence  continues 
under  digital  control  until  the  entire  scan  is  completed. 

3.6  Experimental  Procedure 

To  do  an  experiment,  it  is  first  necessary  to  obtain  an  ultrahigh 
vacuum.  The  entire  system  is  baked  to  300*C  for  approximately  12  hours. 
After  about  24  hours,  the  system  has  cooled  to  room  temperature  and  the 
vacuum  is  typically  5  *  10  ^  torr.  At  this  point,  however,  there  is 
still  considerable  gas  coming  from  the  large  surface  area  of  the  porous 
adsorbent  in  the  cryopump.  Upon  cooling  the  cryopump  to  LN.,  temperature, 
it  ceases  to  be  a  source  and  oeglns  to  pump.  The  pressure  goes  to 
2  x  10  ^  torr.  Once  the  pump  is  cooled  to  liquid  helium  temperature  the 
pressure  drops  below  10  ^  torr.  By  periodically  refilling  the  helium  Dewar 
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Fig.  11.  Functional  block  diagram  of  automated  signal  integrator  and  digitizer 
showing  Integrator  amplifier  (1),  stepping  motor  and  angle  Indicator 
(S),  and  digital  volt  meters  (DVM) . 
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(at  Che  race  of  20  licars  every  4  days)  ic  is  possible  Co  keep  che  syscem 
below  10  ^  corr  for  several  months.  The  syscem  is  now  ready  for  an 
experiment. 

The  surfaces  we  are  inceresced  in  are  formed  by  cleavage  in 

o 

ulcrahigh  vacuum.  The  cleaving  bic  is  a  30  wedge  ground  from  cool 
sceel.  Ic  is  connecced  Co  che  oucside  by  a  flexible  mecal  bellows  coupling. 

A  sharp  cap  of  che  bic  againsc  che  side  of  che  crystal  is  sufficient  to 
produce  a  clean,  relatively  flat  surface  which  gives  reproducible  data 
The  surfaces  obtained  are  quite  immune  to  contamination.  Nonetheless, 
for  che  mosc  pare,  a  new  surface  was  cleaved  every  day.  This  presented  no 
problem  as  che  crystal  could  be  advanced  along  its  length  in  che  holder, 
permitting  as  many  as  15  cleavages  per  crystal. 

To  proceed  with  a  run,  che  pressure  controller  in  che  gas  handling 
syscem  is  activated  to  provide  a  constant  beam  Intensity.  The  intensity 
of  che  beam  can  be  measured  with  che  scattering  detector  by  moving  Che 
sample  ouc  of  the  way.  The  sample  is  then  returned  to  the  center  of  the 
vacuum  system,  and  che  incident  angle  and  azimuth  are  sec.  The  daca 
handling  syscem  can  then  be  activated,  and  the  experiment  runs  automatically, 
scanning  che  region  from  50°  to  180°  from  the  beam.  We  are  constrained 
by  the  geometry  of  che  vacuum  syscem  from  measuring  che  scattered  flux 
at  less  chan  50°  from  the  beam.  When  a  run  is  completed,  in  approximately 
30  minutes,  new  incident  conditions  are  sec  and  another  run  begins.  The 
data  collected  on  paper  tape  are  then  fed  into  a  separate  computer  and  plots 
are  generated. 
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CHAPTER  4 

HELIUM  SCATTER  INC  FROM  MgO  (OOI) 

4.1  Experimental  Results:  Helium  Scattering 

As  a  background  for  our  examination  of  molecular  diffraction,  we  will 
present  data  for  the  scattering  of  He  from  the  (001)  plane  of  MgO.  Both 
the  [110]  and  the  [100]  incident  azimuths  have  been  investigated.  Only 
the  in-plane  scattering  data  are  presented,  as  these  are  the  ones  of 
interest  for  comparison  with  our  subsequent  studies  on  hydrogen.  All  of 
the  data  in  this  section  were  taken  with  both  the  beam  and  the  surface  at 
300  K.  At  these  low  energies,  internal  excitation  of  the  helium  atoms  is 
forbidden.  The  atoms  should  scatter  as  structureless  particles.  This  will 
provide  us  with  a  reference  to  compare  with  the  scattering  of  hydrogenic 
molecules,  where  we  expect  the  internal  degrees  of  freedom  to  be  important. 

4.1.1  fllO]  Azimuth 

Data  for  several  incident  angles  along  the  [100]  azimuth  are  shown 
in  Figs.  12-17.  The  differential  scattering  Intensity  is  plotted  against 
the  angle  between  the  incident  beam  and  the  detector.  As  previously 
mentioned,  only  that  portion  of  the  distribution  scattered  to  T[>  50*  is 
experimentally  accessible.  All  data  in  this  thesis  have  been  normalized 
to  a  direct  beam  Intensity  of  1000  units.  Measurements  were  taken  at 
0.5  degree  intervals.  In  the  figures,  every  data  point  is  marked  by  a  "V, 
with  a  line  connecting  each  point.  To  reduce  the  noise  to  a  low  level, 
all  data  were  collected  using  the  integrator  described  in  Chapter  3. 

It  is  apparent  that  the  in-plane  distribution  consist  of  a  strong 
specular  peak,  flanked  by  broad  diffraction  lobes  on  either  side.  Recall 
that  the  specular  peak  occurs  at  7]  ■  180  -  29^.  Positive  order  diffraction 
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Fig.  12.  H«lium  scattering  from  MgO  (001),  [110]  azimuth.  30.38*  incident  angle 
Specular  intensity  ■  183  units.  Data  normalized  to  1000  units  direct 
beam  intensity.  Data  points  every  .45  degrees,  each  connected  with 
line.  Seam  and  surface  at  300  K.  Expected  peak  position  denoted  by 
arrow. 
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Ftg.  13.  Helium  scattering  from  MgO  (001),  [110]  azimuth,  9 
degrees.  Specular  intensity  -  144  units.  0 
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Fig.  14.  Helium  scattering  from  MgO  (001),  (110]  azimuth.  45.0*  incident 
angle.  Seam  and  surface  at  300  K.  Specular  intensity  "  34.9 
units.  Direct  beam  *  1000  units. 
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Fig.  17.  Helium  scattering  from  MgO,  (001),  [110]  azimuth.  Incident 
angle  “  58.7  degrees.  Specular  intensity  -  44 . 3  units. 
Complete  regular  scan  taken  at  10  sec.  per  point,  10  sec. 
time  constant;  inset,  low  noise  blow  up  of  selected  portion 
showing  the  small  peak  near  7)  •  70*.  Taken  at  60  sec.  per 
point,  60  sec.  time  constant  on  integrator. 
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peaks  will  occur  ac  angles  T)  greater  chan  the  specular  angle;  negative 
order  diffraction  peaks  will  occur  at  T)  less  chan  the  specular  angle. 

At  first  glance  these  lobes  appear  to  be  simply  the  first  order  diffraction 
peaks.  This  is,  however,  incorrect.  The  expected  ositions  of  the  low 
order  diffraction  peaks  are  also  Indicated  in  the  figures.  As  discussed 
in  Chapter  2,  it  is  not  sufficient  to  calculate  the  angle  to  which  the 
most  probable  wavelength  is  scattered.  The  positions  shown  in  the  figures 
were  obtained  from  peak  shapes  calculated  for  a  beam  with  a  velocity 
distribution  conforming  to  Eq.  (7),  assuming  that  the  transition  proba¬ 
bility  is  independent  of  wavelength.  Despite  that,  a  considerable  dis¬ 
crepancy  between  the  positions  calculated  in  this  manner  and  the  peak 
positions  actually  observed  is  apparent.  The  experimental  peaks  in  all 
cases  occur  farther  from  the  specular  direction  than  expected.  In  fact, 
for  several  of  the  curves,  the  predicted  peak  lies  near  the  bottom  of  the 
valley  between  the  specular  and  the  diffraction  peak.  In  addition,  at  high 
d0,  that  is,  toward  normal  incidence,  there  is  evidence  for  a  small  additional 
peak  at  the  bottom  of  the  valley.  Although  this  peak  is  only  slightly  above 
the  noise  level,  it  is  quite  reproducible.  It  can  be  seen  for  all  Incident 
angles  greater  chan  50*;  below  that,  it  disappears  into  the  baseline. 

The  additional  peak  can  be  seen  near  7)  -  70*  in  Figs.  16  and  17. 

To  examine  this  anomalous  peak  more  closely,  several  scans  of  the 

selected  area  were  taken  using  Che  integrator  at  its  maximum  setting  to 

optimize  signal  to  noise.  The  typical  runs  shown  in  this  chapter  were  taken 

using  *»  1  mv  sensitivity  on  the  lock-in,  and  integrating  for  10  seconds  at 

I 

each  angle,  with  a  10  second  time  constant  on  the  integrator.  For  these 
short,  low  noise  scans,  the  lock-in  was  operated  at  50  or  100  uV  sensitivity, 


while  integrating  for  60  seconds  at  each  point,  with  a  60  second  time 
constant  on  the  integrator.  With  this  improvement  in  signal  to  noise, 
the  detailed  structure  in  the  scattered  distribution  can  be  seen  clearly. 

An  example  is  shown  in  the  inset  in  Fig.  17.  The  additional  peak  is 
clearly  not  due  to  background  noise. 

4.1.2  flOO)  Azimuth 

R.  Grant  Rowe  has  previously  presented  data  for  the  scattering  of  He 
from  MgO  (001)  along  the  [100]  incident  azimuth.^  Neither  a  displacement 
of  the  diffraction  peaks  from  the  expected  positions,  nor  any  additional 
peaks  were  observed.  In  an  attempt  to  clarify  the  situation,  measurements 
along  the  [100]  were  repeated  on  the  same  surfaces  on  which  the  [110]  data 
above  were  obtained.  For  the  most  part,  the  [100]  and  [110]  data  were 
taken  in  pairs  for  a  given  Incident  angle,  a  measurement  on  one  azimuth 
immediately  following  the  other  on  the  same  surface. 

Along  the  [100],  the  in-plane  diffraction  peaks  have  indices  (n,n); 
the  lattice  periodicity  along  this  direction  is  4.20  A.  In  Figs.  18-21  the 
data  for  the  [100]  azimuth  can  be  seen  for  comparison  with  the  data  along 
the  [110]  presented  earlier  in  this  chapter.  The  expected  positions  of  the 
first  order  diffraction  peaks  is  also  shown.  For  the  [100]  the  agreement 
between  the  expected  and  observed  positions  is  much  better  than  for  the 
[110].  The  [100]  data  presented  here  also  agree  well  with  previous  data 
on  this  surface.  The  anomalies  observed  with  He  diffraction  along  the  [110] 
azimuth  are  clearly  not  artifacts  of  the  experiment.  Rather  they  must  be 
Indicative  of  the  interactions  between  He  and  MgO. 


There  are  two  aspects  of  the  [110]  data  which  require  examination; 
a)  the  systematic  displacement  of  the  peaks  from  the  expected  positions, 
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Fig.  20.  Helium,  MgO  (001),  (100)  incident  azimuth.  9  »  51.4*. 

65. 5  units  speculsr  intensity.  To  be  compared  to  Fig.  IS 
for  (110)  azimuth. 
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Fig.  21  Helium  scattering  from  MgO  (001),  {100]  azimuth.  8  •  60.5 

degrees.  38.8  units  specular  intensity.  To  be  compared  with 
Fig.  16  for  (110]  azimuth. 
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and  b)  Che  presence  of  small  additional  peaks.  Neither  effect  is  apparent 
on  Che  [100]  azimuth.  The  small  peaks  do  not  appear  at  the  j  order 
diffraction  positions.  This  might  be  expected  if  Che  surface  were  to 
reconstruct  to  a  periodicity  twice  the  bulk  lattice.  The  displacement 
could  result  from  a  straightforward  superposition  of  higher  order  diffrac¬ 
tion  peaks.  However,  since  the  expected  first  order  peak  often  occurs  near 
the  bottom  of  the  valley,  a  very  large  intensity  in  the  higher  orders 
relative  to  the  first  order  would  be  required.  This  appears  unlikely. 

A  possible  explanation  for  the  peak  displacement  is  a  diffraction 
probability  which  is  strongly  a  function  of  wavelength.  Since  the  high 
energy,  short  wavelength  particles  are  diffracted  closest  to  the  specular, 
it  is  necessary  to  severely  attenuate  these  in  order  to  move  the  peak  out¬ 
ward.  If  these  particles  were  preferentially  diffracted  into  the  higher 
order  peaks,  the  diffraction  profile  would  shift  outward.  Furthermore,  if 
the  transition  probability  were  not  monotonic,  but  oscillatory  in  wavelength, 
this  might  produce  the  small  anomalous  peak.  It  remains  to  make  these 
arguments  quantitative.  To  this  end  we  have  carried  out  calculations  of  the 
diffraction  intensity. 

4.2  Semlc lass  leal  Calculations  of  Scattering 

We  wish  to  calculate  diffraction  probabilities  as  a  function  of  wave¬ 
length,  in  such  a  manner  as  to  facilitate  comparison  with  the  experimental 

data.  The  semic lass  leal  approximation  has  proved  useful  in  calculations  of 

18  19 

low  energy  elastic  as  well  as  Inelastic  collisions  in  the  gas  phase,  '  ’  and 
has  recently  been  extended  to  elastic  collisions  at  surfaces.^0 

Basically,  the  approach  is  to  use  classical  mechanics  for  the  collision 
dynamics,  but  quantum  mechanical  superposition  of  the  resulting  scattering 
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Using  Eq.  (16),  this  implies  in  one  dimension  that  | |  «  1,  that  is,  the 
wavelength  does  not  change  appreciably  over  a  distance  comparable  to  the 
wavelength.  In  other  words,  over  small  distances  the  motion  is  plane  wave 
like.  Under  these  conditions,  Eq.  (12)  becomes 


2mT  (VS(r))2 

*  3 


E-V(r) 


(18) 


This  is  just  the  Hamilton- Jacobi  equation  of  classical  mechanics,  where 
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S(r)  is  Hamilton's  characteristic  function.  Equation  (18)  can  be  inte¬ 
grated  to  give 


r 

S(r)  -  p  .  r  +  p  •  ds  ,  (19) 

rl 

the  integral  being  along  the  classical  path  through  the  potential.  We 
see  that,  in  the  classical  limit,  the  "phase"  of  our  wavefunctlon  is  just 
the  action  along  the  classical  path  in  units  of  h.  In  a  straightforward 
manner  we  can  derive  a  relation  for  the  amplitude  of  the  wavefunctlon  using 
the  continuity  equation  Eq.  (16),  to  yield 
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The  term  in  brackets  is  the  Jacobian  of  the  transformation  from  the  final 


coordinates  (x,y)  to  the  initial  coordinates  (x0>yQ)* 
Our  semic lasslcal  wavefunctlon  thus  becomes 


♦  (r) 
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(21) 
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Ic  remains  to  use  quantum  superposition  in  the  proper  manner  to  find  the 
scattering  amplitudes.  In  general  there  will  exist  several  classical 
paths  through  the  potential  connecting  the  initial  and  final  states.  All 
of  the  quantum  effects  in  the  semiclassical  approximation  come  from  the 
superposition  of  these  trajectories. 

We  will  write  the  wavefunction  in  a  Fourier  series  with  coefficients 

Fn  «(,)  ■TTjj’  *  (x,y,*)  exp(-i(K  +  2  >-t)dxdy  ,  (22) 

n  ,m  a  •  *  *  on, in 

x  y  unit  cell 

where  we  have  made  use  of  the  periodic  nature  of  the  surface  potential 

The  diffraction  intensities  will  be  given  by  v  | F  j*.  Substituting  our 

z  n,m 

semiclassical  wavefunction,  we  obtain 
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(23) 


We  are  interested  in  Eq.  (23)  only  in  the  limit  z-4®,  that  is,  far  from 

Che  surface.  The  integral  in  Eq.  (23)  is  over  (x,y),  that  is,  over 

different  starting  points  for  the  trajectories.  It  is  thus  an  integral  over 

all  the  possible  trajectories.  If  S/&  is  large,  the  phase  of  the  integrand 

varies  rapidly  as  a  function  of  (x,y).  Integrals  of  this  type  are  usually 
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evaluated  in  the  stationary  phase  approximation.  In  one  dimension,  the 
method  of  stationary  phase  applies  to  integals  of  the  form 

_  f  ,  .  if(x). 

I  •  J  g(x)«  dx 


where  f(x)  is  a  rapidly  varying  function  of  x. 
approximation  is 


)* 


If  (x  t) 


The  stationary  phase 


(24) 
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where  the  x^  ere  solutions  to  •  0,  that  Is,  they  are  the  points  of 

stationary  phase.  The  stationary  phase  condition  applied  to  Eq.  (23), 
however,  is  precisely  the  diffraction  condition,  Eq.  (3), 


K  •  K  +  G 


n,m 


(3) 


The  diffraction  condition  is  thus  obtained  in  a  very  natural  way  within 
the  semiclassical  approximation. 

We  have  concluded  that  only  those  points  on  the  surface  which  lie  along 
the  classical  trajectories  will  contribute  to  scattering,  and,  furthermore, 
that  only  those  trajectories  for  which  Eq.  (3)  is  satisfied  will  contribute. 
The  primitive  semiclassical  approximation  to  the  scattering  Intensities 
results  from  a  straight  quantum  superposition  of  the  several  (two)  contri¬ 
buting  amplitudes.  We  thus  have  the  primitive  semiclassical  approximation 
for  the  scattering  amplitude 


where 


F 

n  ,m 


-  tr/u.: . 

3(p  ,p ) 


(25a) 


9  ■  -  x£xdt  ♦  ypydc  +  *P2dc  (25b ) 

and  the  sum  is  over  all  the  classically  allowed  trajectories  in  the  unit 
cell  satisfying  Eq.  (3),  the  diffraction  condition. 

The  stationary  phase  formula,  Eq.  (24),  was  derived  by  expanding  the 
phase,  f,  in  a  quadratic  expression.  It  becomes  invalid  if  f"  is  small. 

This  will  occur  whenever  two  points  of  stationary  phase,  x^,  and  x^ ,  coalesce. 
The  trajectories  for  which  this  occurs  are  of  considerable  interest  in 
scattering  calculations.  To  avoid  the  singularity  which  occurs  in  using 
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Eq.  (24)  when  f"  is  small,  the  uniform  semiclassical  approximation  was 

developed. By  expanding  f(x)  near  the  points  of  stationary  phase  as 

a  cubic,  a  new  stationary  phase  formula,  involving  Airy  functions,  can  be 

derived  which  is  well  behaved  in  this  critical  region.  The  scattering 

22  19b 

intensities  in  the  uniform  semiclassical  approximation  become  * 

Pn,m  ‘  <KlMn2l*)2  nq*Al2(-q)  +  (|ni|*-jn2|i)TTq2Bi2(-q)  (2 

where  q  -  (3/41©^  -  ]  ) 2 ^ 3  (2 


n‘ '  TljvfZ 
a<w 

The  phases  cp^  for  the  two  trajectories  are  given  by  Eq.  (25b).  We  will  use 

the  semiclassical  approximation  in  this  form. 

There  is  an  alternative  derivation  of  the  semiclassical  approximation 

which  proceeds  not  from  the  Schrodinger  equation,  but  from  the  Feynman 
43 

path  integral  formulation  of  quantum  mechanics.  It  is  mentioned  here 
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contribute  to  the  scattering  amplitude,  each  with  its  own  phase.  All 
that  is  necessary  to  get  to  the  seroic lassical  approximation  is  to  evaluate 
this  integral  in  the  stationary  phase  approximation;  now  only  the 
classically  allowed  paths  contribute.  All  the  approximations  made  in 
the  earlier  derivation  are  thus  equivalent  to  only  one  common  assumption, 
that  is,  the  validity  of  the  stationary  phase  approximation. 

4.2.2  Calculation  of  Diffraction  Intensities 

We  wish  to  calculate  the  diffraction  probabilities  using  Eqs.  (26)- 
(28).  We  desire  a  potential  which  contains  the  essential  features  of  the 
surface  potential  yet  which  is  caiculationally  simple.  Only  in  plane 
scattering  from  a  potential  with  one  dimensional  periodicity  will  be 
considered.  The  potential 


V (x,z)  - 


0 


z  >  & 

i  >  z  >  Zo 
*<Zo 


.2rrx, 

c°s(~> 

,2rut. 

co»(-r> 


(30) 


shown  schematically  in  Fig.  22  has  previously  been  used  by  Steele  and 

22 

Doll  “  for  calculation  of  gas-surface  interactions.  It  is  a  corrugated 

hard  wall  with  a  square  well  in  fronc  of  it.  The  real  surface  potential 

has  a  gentle  long  range  attractive  part  and  a  short  range,  fairly  hard 

repulsion.  The  long  range  attraction,  however,  has  been  hown  to  have 

44  45 

little  effect  on  the  scattering.  ’  Only  the  hard,  short  range 
potential  is  important.  Thus,  although  this  model  potential  is  crude, 
it  contains  the  essential  features  of  the  interaction,  and  has  been 
adopted  for  our  calculations. 
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Our  semiclassical  approximation  should  be  most  appropriate  for  short 

20d 

wavelengths  and  for  gradually  changing  potentials.  Masel,  however, 

has  compared  the  semiclassical  solution  to  an  exact  quantum  solution  for 

normal  incidence  in  a  simple  corragated  hard  wall  potential.  This  hard 

wall  should  be  an  extreme  test  of  the  semiclassical  theory.  Essentially 

exact  agreement  between  the  quantum  and  semiclassical  results  was  obtained 

down  to  wavevectors  lea  -  10,  that  is,  up  to  wavelengths  X  1.8  A. 

Since  this  is  more  chan  three  times  Che  most  probable  wavelength  in  a 

thermal  beam  at  300°K,  the  semiclassical  approximation  should  give 

excellent  results,  even  for  this  simple  potential. 

The  classical  trajectories  in  this  piecewise  constant  potential  are 

straight  lines.  One  such  trajectory  is  shown  in  Fig.  22.  There  are  two 

equivalent  trajectories  per  unit  cell,  corresponding  to  the  two  points  of 

equal  slope.  These  are  denoted  by  A  and  B  in  the  figure.  For  this  case, 

2  2  46 

considerable  simplification  of  Eq.  (26)  results.  We  obtain  ’ 

Pn(9,E)  -  4rq*|Ai(-q)|^  ,  (31) 


where 


n 


“i 


^2  " 


*  .  2 
(1-tan  y  tan  <y)(l+tan  cr) 

•A... 

o  k  a  n 


(32) 


The  angles  are  given  by  Y  ■  90  -  8  and  sin  Y  *  yi+V/E  sin  Y* ,  the 
referring  to  the  corresponding  angle  measured  inside  Che  well,  x^  is 
the  x  coordinate  at  which  the  trajectory  intercepts  the  hard  wall  and  a  is 
Che  slope  of  the  corrugation  at  chat  point.  P^O.E)  is  the  probability  of 
a  particle  of  energy  E  being  diffracted  into  diffraction  peak  n  at 


- 
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Using  the  semiclassical  formalism  just  described  we  have  undertaken 

the  calculation  of  the  diffraction  peak  profiles,  closely  following  the 
22 

procedure  of  Doll.  The  semiclassical  intensities  q(0,E)  for 
-5  <  n  <  5  were  thermally  averaged  over  the  energy  distribution  of  a 
thermal  beam 

F(E)dE  -  (kT)'2E  exp(-E/kT)dE  .  (33) 

20c 

The  intensity  at  any  angle  (except  the  specular)  is  given  by 


1(9) 


2  sin  9 
COS  0  -  cos 


|£  E  F(En)P<9.En) 
o  n 


(34) 


where  E  “  E(n,9  ,9)  is  the  energy  of  the  atoms  diffracted  from  9  to  9 
no  o 

in  diffraction  channel  n,  subject  to  Eq.  (3)  and  conservation  of  energy. 

The  P(9,E  )'s  are  the  uniform  semiclassical  intensities  given  by  Eq.  (31). 
n 

dE 

The  term  in  front  of  the  summation  arises  from  —  for  elastic  scattering. 

Within  this  model  the  calculations  are  straightforward  and  can  be  done  for 

a  range  of  parameters.  The  variable  parameters  in  the  model  are  the 

depth  of  the  square  well,  Vq,  and  the  amplitude  of  the  hard  wall,  Zq. 

Calculated  distributions  for  a  range  of  parameters  for  scattering 

along  the  [110]  are  shown  in  Figs.  23-32,  where  is  given  in  degrees 

Kelvin  and  ZQ  in  units  of  the  lattice  parameter,  a. 

The  effect  of  variation  of  V  at  constant  Z  is  shown  in  Figs.  23-26. 

oo 

The  major  effect  of  is  to  refract  the  trajectories  as  they  exit  the 
well.  This  shifts  the  sharp  cutoff  at  high  and  low  angles  away  from  the 
specular.  These  cutoffs  occur  at  the  so-called  rainbow  angles,  a  carryover 
from  classical  scattering.  Classically,  there  exists  a  maximum  and  minimum 
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Fig.  24.  Semic lassical  scattering  distribution  for  •  59*, 

azimuth.  Helium  from  MgO  (001).  V  *  50  K,  Z  •  0. 

o  o 

Beam  •  300  K. 
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Fig.  25.  Seme  as  Fig.  23  except  Vq  •  100  K. 
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angle  at  which  scattering  from  a  periodic  potential  can  occur;  for  our 
case  this  corresponds  to  trajectories  which  strike  the  inflection  points 
on  the  hard  wall,  the  points  of  maximum  slope.  Trajectories  which 
require  a  deflection  larger  than  this  cannot  occur.  This  classical  factor 
enters  our  semlclasslcal  formalism  through  the  Jacobian,  Cl.  The  sharp 
cutoff  is  entirely  the  result  of  our  hard  wall  potential.  For  a 
realistic  potential,  the  inclusion  of  complex  (tunneling)  trajectories 
softens  this  sharp  edge.  For  our  hard  wall  potential,  no  such  complex 
trajectories  exist.  Rainbow  scattering,  however,  is  not  just  a  classical 
effect.  The  same  cutoffs,  albeit  not  as  sharp,  appear  in  a  full  quantum 
treatment  of  surface  scattering. 

We  are  for  the  most  part  interested  only  in  that  portion  of  the 

scattering  distribution  closer  to  the  specular.  As  Vq  appears  to  have 

little  effect  on  this  part  we  will  continue  these  calculations  with 

fixed  at  100  K.  This  corresponds  to  a  well  depth  of  8.62  meV,  or 

about  16  per  cent  of  the  particle's  incident  kinetic  energy.  By 

comparison,  the  Morse  potential  well  depth  for  helium  on  L1F ,  derived 
49 

from  experiment,  is  7.59  meV.  The  MgO  well  depth  can  be  expected 

to  be  of  the  same  order  of  magnitude. 

Variation  of  Z Q  has  a  significant  effect  upon  the  calculated 

profiles.  Shown  In  Figs.  27*30,  as  well  as  Fig.  25,  are  the  distributions 

calculated  for  various  2  for  9  •  59°  and  V  -  100  K,  as  well  as  the 

o  o  o 

experimentally  observed  positions  from  Fig.  17.  It  is  spparent  that  the 
broad  lobes  observed  in  Che  experiments  are  the  sum  of  up  to  5  dif¬ 
fraction  peaks.  Significant  shifts  in  the  peaks  take  place  as  a  function 
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Fig.  27.  Semic lass  leal  scattering  calculations  for  helium,  Z  ■  .05  A, 

o 

Vq  -  100  K,  •  59*.  For  comparison  with  experimental  data  in 
Fig.  17.  The  specular  position  indicated  by  dashed  line. 
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Fig.  28.  Semic lsssical  scattering  distribution*  for  helium,  Z  •  .1  A, 

o 

■  100  K,  •  59*.  Lattice  spacing  •  2.97  A.  Average  over 

thermal  beam  at  300  K,  diffraction  channels  |n|<5.  Experimental 
peak  positions  taken  from  Fig.  17. 
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Scattering  Angle,  rj  (degrees) 

Fig.  29.  Calculated  semic lassical  distributions  forfeiture,  (110]  azimuth 
MgO  (001).  9^  •  59*,  ■  100  K,  ■  .2  A.  Thermal  averaged 

profile  for  Ig€4B  •  300  K.  Compare  to  data  in  Fig.  17.  Experiment 

peaks  indicated  by  arrows. 
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Fig.  30.  Calculated  seralclaf slcal  diffraction  profile  at  59*  Incident 
angle.  -  0.16  A,  •  100  K.  Thermal  beam  at  300  K. 

Experimental  peaks  Indicated  by  arrows.  See  Fig.  17. 

Only  portion  "]  )  50*  Is  experimentally  accessible,  arbitrary 
units.  (110]  azimuth. 
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of  Z  ;  for  large  Z  ,  Che  lobe  spiles  Into  several  peaks.  Significant 
o  o 

oscillations  also  appear  in  the  intensity  of  the  individual  diffraction 
channels.  This  can  be  seen  in  Fig.  31  where  the  diffraction  lobes  have 
been  broken  down  into  their  component  peaks.  These  oscillations  are  due 
to  oscillations  in  the  Airy  function  factor  in  Eq.  31.  Physically,  this 
results  from  the  interference  of  the  two  classical  trajectories  which 
contribute  in  the  semiclassical  approximation  to  a  given  diffraction 
channel. 

A  few  notes  of  caution  on  these  calculations  should  be  made.  We 
have  assumed  chat  there  is  no  classical  multiple  scattering;  Che 
particle  hits  Che  wall  only  once.  This  will  not  be  a  good  approximation 
at  glancing  incidence  angles,  glancing  exit  angles,  or  at  high  Zq.  For 
our  potential,  however,  refraction  at  the  square-well  keeps  the 
trajectories  Inside  the  well  closer  to  normal  incidence.  Furthermore, 
it  is  Che  long  wavelength,  low  energy  particles  for  which  these  glancing 
exit  angles  occur.  These  are  precisely  those  for  which  the  semiclassical 
approximation  is  least  valid.  The  section  of  the  distributions  far  from 
the  specular  is  thus  the  least  accurate.  Fortunately,  it  is  the  phenomena 
which  occur  far  from  this  region,  near  the  specular,  we  are  interested  in. 

4.3  Discussion 

Qualitatively,  both  the  anomalous  features  in  the  [  110]  data  can  be 

explained  at  3  *  59°  with  V  ■  100  K  and  Z  *0.16  k.  The  small  peaks 

o  o  o 

near  *'  •  55°  and  Tj  »  70°  generally  disappear  below  Zg  •  0.14  k.  Similar 
calculations  at  ?o  *  30°  and  Sg  •  51°  for  these  potential  parameters  are 


Calculated  semic last  leal  diffraction  peak  profiles  for  MgO  (001) 

[110]  azimuth.  Helium  scattering  averaged  over  thermal  beam  at 

300  K.  The  individual  diffraction  peaks  up  to  | nj  >3  are  shown 

for  V  -  100  K,  Z  »  .2  A,  at  9  -  59*.  The  total  profile,  as 
oo  o 

displayed  in  previous  figures,  is  the  sum  of  these  individual 
distributions. 
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shown  in  Figs.  32  and  33,  respectively.  These  are  to  be  compared  to  the 
data  in  Figs.  12  and  15.  In  Fig.  32  for  8q  ■  30°,  the  calculated  and 
exrerimental  positions  agree  well;  the  peak  is  somewhat  broader  than 
experimentally  observed.  The  shape  of  this  peak  is,  however,  quite  sen¬ 
sitive  to  V^.  At  •  51°  in  Fig.  33,  the  large  peak  at  1)  ■  99°  fits 
well  to  the  data;  the  small  peak  near  11  ■  89°  has  disappeared  into  the 
inelastic  background.  The  low  angle  lobe  does  not  fit  nearly  as  well, 
however.  The  calculated  peak  occurs  at  T,  •  52°,  7°  below  the  experimental 
peak  at  H  •  59°.  In  contrast,  the  position  of  the  (1,0)  peak  calculated 

using  a  constant  transition  probability  as  described  earlier  gives  a 
o  o 

peak  at  1]  •  63  ,  4  to  the  other  side  of  the  experimental  peak. 

Unfortunately,  we  cannot  study  these  peaks  over  a  wider  range  of  6  ; 

o 

above  8^  *  50°,  the  (1,0)  peak  is  experimentally  inaccessible. 

The  [100]  data  from  Figs.  18-21  do  not  show  any  serious  shifts  of 
extra  peaks.  In  these  semiclass  leal  calculations,  the  scattering  profile 
is  a  sensitive  function  of  the  corrugation  amplitude  T.q  and  the  lattice 
spacing  a;  the  well  depth  has  a  lesser  effect.  Along  the  [ 100]  the 
effective  lattice  spacing  is  4.2  k,  that  is,  v^2  times  longer  than  the 
[llO]  spacing.  If  we  take  Zq  to  be  the  same  along  both  azimuths,  we  can 
qualitatively  consider  the  effect  of  changing  the  lattice  dimension.  As 
£  Increases,  or  ,  the  maximum  slope  of  the  potential  wall,  decreases. 

The  rainbow  angles,  8_  ,  move  inward  according  to  ■  8  +  a  .  A 

major  effect  will  thus  be  to  concentrate  the  scattering  into  a  narrower 
angular  region;  all  other  things  equal,  we  would  thus  expect  the 
scattered  distribution  to  become  sharper  and  taller. 
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Fig.  32.  Calculated  semic lassica 


asslcal  distribution  for  helium  scattering 
•  .16  A  and  V  -  100  K.  The  experimental 


at  9  -  30*.  Z 

o  o 

peak  position  from  Fig.  12  Is  indicated.  The  small  bump  near 
■  140*  Is  the  (1,0)  diffraction  peak  which  is  energeticly 
closed  to  most  of  the  beam  molecules. 
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Fig.  33.  Semic lassical  distribution  for  3  -  SI*,  Z  *  .16  A,  V  -  100  K. 

o  0 

For  comparison  with  experimental  data  In  Fig.  15.  Experimentally 
observed  peaks  indicated  by  an  arrow.  Only  the  portion  at 
50*  is  experimentally  accessible. 
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Semiclass  leal  calculations  have  been  done  along  Che  r100]  azimuth 
using  Che  same  procedure  as  for  the  [110].  It  was  not  possible,  however, 
to  fit  Che  experimental  peak  positions  to  the  calculations  using  the  same 
potential  parameters  as  determined  for  the  [110].  The  same  corrugation 
amplitude,  ■  0.16  A,  was  used  but  a  much  deeper  well,  Vq  »  200  K  was 
required.  Calculated  distributions  using  these  parameters  for  ■  60° 
and  ■  30°  are  shown  in  Figs.  34  and  35  for  comparison  wich  the  data 
in  Figs.  21  and  18.  The  positions  of  the  major  peaks  fit  well.  However, 
for  ■  60°  in  Fig.  34,  the  small  peak  near  7|  ■  67°  is  not  visible  in 
Che  data  even  though  a  similar  peak  was  clearly  visible  in  the  [110] 
data.  We  can  only  conclude  that  it  has  become  lost  in  the  inelastic 
tail  of  the  specular  peak.  The  peak  at  T|  •  96°  for  0Q  •  30°  in  Fig.  35 
also  fits  well  to  the  data.  The  small  lobe  to  the  right  of  the  specular 
angle,  near  140°,  is  the  high  energy  tail  of  the  (1,0)  diffraction  peak. 
This  channel  is  energetically  closed  to  a  significant  fraction  of  the 
beam  molecules.  This  small  lobe  is  not  observed  in  the  data. 

It  should  not  be  surprising  that  a  different  potential  is  required 
along  Che  [100]  chan  along  Che  [110].  As  shown  in  Fig.  1,  the  [110]  is 
the  direction  between  nearest  like  ions  on  the  surface,  whereas  the  [100] 
direction  passes  through  alternate  positive  and  negative  ion  cores. 
Tsuchida^  has  calculated  the  helium-LiF  potential  by  summation  of 
pairwise  interactions.  The  He-MgO  interaction  should  be  qualitatively 
similar.  The  V  -  +  1000  cal/mole  equipotential  surface  of  this  potential 
is  shown  schematically  in  Fig.  36.  One  quarter  of  the  unit  cell  is  shown. 
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Fig.  34.  Semiclasaical  calculation  for  [100]  azimuth.  9  ■  60* 

•  ® 

Z  *  0.16  A,  V  •  200  K.  Experimental  peak  denoted 
o  o 

by  arrow. 
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Fig.  35.  Semic lass  lea l  calculation  for  fiOOJ  axlauth.  Experimental 

peak  denoted  by  arrow.  Z  -  0.16  A,  V  -  200  K,  9  -  30* 

o  o  o 
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Schematic  equi. potent  ial  surface  (V  •  +  1000  cal/mola)  for 
helium  on  LLF  (001),  as  calculated  by  Tsuchlda  (ref.  SO) 
Contours  denote  distance  from  surface  In  angstrom a. 
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The  contours  represent  the  distance  from  the  surface  plane  to  the 
equipotential.  Since  the  mean  energy  in  a  beam  at  300  K  is  1192  cal/mole, 
this  surface  approximates  the  classical  turning  point  surface  and  can  be 
roughly  compared  to  our  hard  wall.  Along  the  [ 110]  the  equipotential 
appears  to  be  a  simple  sinusoidal  function  of  distance,  with  the  valley 
midway  between  the  like  ions.  In  fact,  the  effective  roughness  of 
Tsuchlda's  potential,  similar  to  our  Z^,  is  0.17  A  along  the  [110*,  in 
excellent  agreement  with  our  0.16  A.  This  is,  however,  perhaps  fortuitous. 
Along  the  [100]  azimuth,  the  potential  is  more  complex  with  2  maxima  per 
unit  cell.  Thus  it  may  be  that  our  model  potential  is  too  simple  to 
handle  scattering  along  the  [100]  azimuth.  A  more  realistic  potential, 
such  as  Tsuchlda's,  could  be  handled  by  the  same  semic lassical 
formalism  we  have  used.  It  would  require,  however,  numerical  integration 
of  the  classical  equations  of  motion  through  Che  potential.  Note  that 
since  there  are  4  inflection  points  per  unit  cell  along  the  ,  110'1  for 
this  potential,  there  will  in  general  be  4  contributing  trajectories 
Instead  of  the  two  in  our  calculation.  This  could  significantly  change 
the  Interference  structure  in  the  scattered  distributions.  Although 
such  a  calculation  could  be  done,  it  would  be  extensive  and  it  thus 
would  defeat  one  of  our  main  purposes,  that  is,  to  have  as  simple 
calculatlonal  scheme  to  aid  in  experimental  interpretation. 

The  potential  parameters  determined  here  should  be  considered  only 
an  approximation.  The  model  is  not  Intended  to  be  quantitative.  Although 
the  calculation  works  well  for  most  of  the  data,  the  position  of  the 
low  7)  peak  at  •  51°  on  the  [110]  as  well  as  the  absence  of  the  small 
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peaks  in  Che  daca  9^  ■  60°  on  Che  [100^  remain  troublesome .  This  should 

noc,  however,  be  due  co  a  deficiency  in  Che  semiclassical  approximation 

buC  racher  due  co  inadequacies  in  che  model  pocencial  and  the  artificial 

limitation  to  in-plane  scattering.  The  primary  purpose  of  these 

calculations  has,  however,  been  served;  thac  is  co  qualitatively  explain 

Che  experimental  peak  positions  and  fine  scructure.  We  conclude  Chat 

Che  transition  probability  for  elastic  diffraction,  T  (X),  is  a  strong 

function  of  wavelength.  In  addition,  there  is  considerable  intensity 

in  the  higher  order  diffraction  peaks  up  to  |n|  -  5.  It  is  the 

superposition  of  these  peaks,  each  having  its  own  structure  dictated  by 

T  (X),  which  results  in  the  experimentally  observed  distributions, 
n 

A  comment  on  the  use  of  a  thermal  beam  for  elastic  diffraction 
studies  is  in  order.  Our  calculations  indicate  significant  diffraction 
intensity  out  to  |n  ■  5.  Noc  only  are  the  individual  diffraction  peaks 
broad  due  to  the  distribution  of  velocities  in  the  beam  distribution, 
but  there  is  also  extensive  peak  overlap.  The  fine  structure  and 
intensities  of  che  individual  peaks  are  thus  not  observable  with  a  thermal 
beam.  Considerable  information  is  lost.  It  is  clear  that  quantitative 
studies  of  the  gas-surface  potential  through  elastic  diffraction  of 
molecular  beams  require  a  monoenergetlc  beam. 


CHAPTER  5 


SCATTERING  OF  HYDROGEN IC  MOLECULES 
5 . I  Experimental  Results:  HD 

In  this  chapter  we  will  present  the  scattering  of  the  hydrogenic 
molecules  H, ,  D, ,  and  HD  from  the  (001)  plane  of  MgO  along  the  (110J 
incident  azimuth.  We  wish  to  investigate  the  participation  of  rotational 
modes  of  simple  molecules  in  gas  surface  scattering.  Particular  emphasis 
will  thus  be  placed  on  the  scattering  of  HD,  as  the  rotational  modes  of 
this  molecule  most  strongly  couple  to  the  surface.  Data  were  taken  with 
the  surface  at  300  K.  for  beam  temperature  at  both  300  K  and  80  K. 

The  results  of  scattering  experiments  with  HD  at  a  beam  temperature 
of  300  K  are  shorn  in  Figs.  37*41.  Upon  comparison  with  the  helium 
measurements  presented  earlier,  several  differences  are  immediately 
apparent.  There  is  considerably  more  inelastic  scattering  for  HD  than 
for  the  helium;  the  specular  peak  has  been  severely  attenuated,  and  the 
diffraction  lobes  are  less  well  defined.  However,  despite  the  large 
background  due  to  thermal  inelastic  events  and  peak  overlap,  the  features 
are  clearly  marked.  Instead  of  the  smooth,  broad  diffraction  peaks  which 
we  observed  for  helium,  the  HD  distribution  is  dominated  by  a  series  of 
sharp  and  very  narrow  peaks.  It  is  immediately  clear  that  these  peaks 
result  from  a  different  process  than  the  essentially  elastic  diffraction 
observed  for  helium  scattering.  Recalling  our  discussion  in  Chapter  2, 
it  appears  these  peaks  are  due  to  rotationally  inelastic  diffraction 
processes. 

Of  particular  Interest  in  these  distribution  is  the  sharpness  of  many 
of  the  peaks.  Although  only  the  tops  are  visible,  they  appear  to  be 
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HD  scattering,  MgO  (001),  (110)  azimuth,  300  K  be 
9.2*  incident  angle.  Specular  -  26.3  units. 


Fig.  40.  HD  scattering,  9  -  49. 5* ,  (IlOj  azimuth.  Specular  peak  ■ 

11.9  units . 
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considerably  narrower  chan  che  elastic  peaks  we  saw  for  helium.  This 
raises  a  quescion:  why  should  Che  inelastic  peaks  be  narrower  than  the 
elastic  peaks.  This  will  be  explored  later  in  this  chapter.  For  clarity 
and  future  reference,  the  expected  positions  of  the  rotational  loss  peaks 
have  been  labeled  in  these  figures.  The  calculation  of  the  peak  positions 
is  not  straightforward,  however.  It  is  intimately  related  to  the  calcu¬ 
lation  of  the  peak  shapes.  The  labeling  of  the  peak  positions  has  not  been 
justified  at  this  point.  This  entire  matter,  however,  will  be  covered 
extensively  later  in  this  section. 

In  order  to  more  clearly  see  che  shapes  and  positions  of  these 

rotational  peaks,  several  low  noise  scans  have  been  taken  of  selected  areas, 

using  Che  integrator  in  the  manner  described  in  che  previous  chapters. 

Ctte  such  scan  is  shown  in  Fig.  41,  where  the  peaks  for  9  -  13.2*  are 

o 

clearly  revealed.  There  appears  to  be  a  distinct  asymnetry  in  the  peak 
labeled  (l,0)-l  near  71  *  130*.  It  rises  almost  vertically  on  the  low  angle 
side,  with  a  more  sloping  decrease  on  the  right  side. 

The  large  background  in  the  data  presented  in  this  chapter  has  two 
origins;  a)  thermal  Inelastic  events,  which  broaden  the  peaks,  and  b)  peak 
overlap.  Significant  reduction  of  the  thermal  background,  requires  cooling 
Che  crystal  Itself  below  che  ambient  temperature.  This  was  not  possible 
in  these  experiments.  Measurements  have,  however,  been  taken  with  an 
80  K  beam  of  HO,  again  incident  along  the  [110]  azimuth.  These  data  are 
shown  in  Figs.  42-44.  Recalling  that  X  -  h/p«» h/v 2M^kT,  and 

cos6  •  cos9  4*  nX/d  ,  (4) 

o 

we  expect  the  elastic  diffraction  peaks  to  move  outward,  away  from  che 
specular,  as  the  temperature  of  the  beam  is  lowered.  By  lowering  T, 


Scattering  Angle,  ^(Degrees) 


Fig.  6i.  Slow,  low  noise  scsn  of  selected  portion  of  distribution 
showing  shsrp  rotational  peaks.  0^  •  13.2*  for  HD,  [110] 

azimuth.  Data  taken  at  60  sec.  time  constant,  60  sec. 
per  point. 
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Fig.  44.  HD  scattering  for  30  K  beam.  Incident  angle  ■  29.7*. 
[HO]  azimuth. 
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we  can  hope  Co  decrease  peak  overlap  and  more  clearly  resolve  the 
individual  peaks. 

Lowering  che  temperature  of  the  source  has  effects  ocher  than 
changing  X  and  shifting  che  peaks  outward.  Cooling  che  beam  will 
also  shift  the  distribution  of  rotational  states  toward  lower  J,  as 
shown  in  Fig.  45.  Furthermore,  as  T  is  lowered,  chat  is,  as  the  energy 
of  che  beam  is  lowered,  fewer  elastic  scattering  channels  are  open. 

One  effect  of  this  can  be  seen  in  comparing  Fig.  37  to  Fig.  42,  for 
scattering  at  approximately  ^  ■  50*  at  300  K  and  80  K  respectively. 

The  specular  intensity  (l.e.  che  clastic,  non-dif f racted  intensity)  is 
approximately  twice  as  large  at  80  K  as  at  300  K. 

5.2  Extremal  Effects 

Of  particular  interest  in  Che  data  is  che  apparent  sharpness  of  che 
peaks,  as  well  as  their  peculiar  shape.  Furthermore,  it  can  be  noted 
that  although  some  of  che  transitions  are  strong,  others  are  not  observed 
at  all,  for  example,  at  9^  ■  30*,  the  (1,0) +2  is  clearly  visible  but 
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Fig.  45.  Rotational  state  population  for  normal  HD  at  300  K  and 
80  K,  from  Table  l. 


where  again  k0  *  and  k  are  the  incident  and  scattered  wave 

vectors,  respectively,  and  &E  is  the  energy  transferred  from  rotation 
to  translation.  Note  that  for  inelastic  scattering,  the  scattered  wave¬ 
length,  \,  is  no  longer  equal  to  the  Incident  wavelength,  \q.  We  will 
at  present  confine  ourselves  to  in-plane  scattering  along  the  UlO) 
azimuth;  for  this  case  m«0  and  the  Inelastic  diffraction  peaks  are  labeled 
(n,0)  +  l.  From  Eq.  (5)  we  previously  obtained  Eq.  (6)  for  the  scattered 


angles 


cos$  ■  —  (cosS 


5^-  -  ( l  ♦  2-EM  X2/h2f* 
o 


Although  at  first  these  relations  appear  normal,  on  closer  examination, 

peculiar  behavior  is  observed.  Consider  the  case  of  a  loss  transition 

from  the  (1,0)  diffraction  peak,  for  which  both  ^E  and  n  are  negative. 

In  this  instance,  \/\  is  an  increasing  function  of  wavelength,  whereas 
nX 

cos3  ——  is  a  decreasing  function.  Cos§  is  the  product  of  these  two 
o  a 

terms  and  thus  the  relationship  can  be  double  valued,  that  is,  it  may 

have  a  turning  point.  Such  a  relation  between  Vq  and  ~  is  shown  in  Fig.  46, 

where  it  is  important  to  recall  that  6  and  "  are  related  by  T)  •  180-6-6^. 

The  effect  of  this  extremum  on  the  rotational  inelastic  peaks  can  now 

dF 

be  clearly  seen.  We  know  that  the  differential  scattering  intensity  jr 
can  be  written  as 


2<I>I(>.0>«V  35* 


where  KX^)  and  t 0Q)  ere  the  transition  probability  and  the  wavelength 
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Fig. 


In-plane  scattering  of  HD  of  MgO  (001)  azimuth  9  ■  13.2*. 

0 

a)  Scattering  angle  as  a  function  of  \  for  elastic  (1,0) 

o 

peak.  X  is  Che  most  probable  wavelength  in  the  beam, 
mp 

b)  Scattering  angle  as  a  function  of  \  for  the  (1,0)  -  1 

^  MM 

rotational  loss  peak,  showing  the  extremum  at  (\  -  ). 

c)  Extremal  scattering  intensity  for  (!,0)-l  transition, 

M 

with  singularity  at  the  extremal  angle  *  . 
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distribution  respectively.  Z(I)  is  the  population  of  molecules  in  the 
proper  rotational  state  for  transition  I.  It  is  apparent  from  Fig.  46b 

J  p 

chat  -jT-  become  infinite  at  the  turning  point,  that  is,  at  the  angle 

9  •  9  .  This  is  shown  in  Fig.  46c. 
m 

This  simple  argument  predicts  a  singularity  in  the  scattering 
intensity  for  these  rotational  processes  at  certain  scattered  angles,  9 

m 

called  the  extremal  angles.  Ue  will  denote  this  process  as  "extremal 

scattering",  (^lalitative ly  we  can  understand  this  process  by  referring 

again  to  Fig.  46.  At  the  extremal  angle,  a  large  range  of  wavelengths  is 

scattered  to  the  same  scattered  angle.  The  kinematics  of  the  inelastic 

scattering  process  have  acted  to  focus  Che  scattered  particles.  This  is 

in  contrast  to  the  elastic  diffraction  peaks,  where  each  wavelength  goes 

d\ 

to  a  different  angle.  Extremal  scattering  will  occur  whenever  ■  0. 

If  follows  from  Eq.  (6)  that  Che  extremal  angle  is  given  by 


cos  9 

a 


(1 


2.2  4 

n_h__^5 


2d2££M  cos29 
8  ° 


(36) 


The  wavelength  of  the  incident  molecules  scattered  at  8  is 

m 

2 

X  •  nh  /  2  &EM  dcos9  .  07) 

9  gO  '  ' 

Extremal  scattering  peaks  can  thus  only  occur  when  n,  the  diffraction 
index,  and  &E  the  energy  change,  have  the  same  sign  (nonzero).  The  extremal 
scattering  peaks  will  therefore  appear  between  the  elastic  diffraction  peak 
and  Che  specular  (0,0)  peak.  Although  these  considerations  predict  a 
singularity  in  the  differential  scattering  Intensity,  this  singularity 
will  be  removed  by  integration  over  the  finite  aperture  of  the  detector  in 


any  real  experiment.  However,  we  expect  these  peaks  to  remain  quite 
distorted  and  sharp  when  compared  with  the  elastic  scattering  peaks. 

To  get  a  clear  insight  into  the  magnitude  of  extremal  scattering, 
we  have  undertaken  some  calculations  to  simulate  the  experimental  peak 
shapes.  Calculations  were  done  for  two  different  beam  temperatures, 

300  K  and  30  K.  The  elastic  diffraction  peaks  will  shift  outward,  away 
from  the  specular,  as  the  temperature  is  decreased  and  the  wavelengths 
increase.  The  positions  of  the  extremal  peaks,  however,  given  by  Eq.  (36), 
depend  only  on  and  <1E  and  should  not  change.  The  shapes  should  change 
however,  as  the  wavelength  distribution  shifts  through  the  singularity. 

Several  factors  should  influence  the  observed  shapes  and  intensities 
of  the  rotational  peaks.  Briefly,  the  most  important  are: 

1)  The  transition  probability,  T(X),  as  a  function  of  wavelength. 

2)  The  initial  wavelength  distribution,  f(X),  in  this  case  a  modified 
Maxwe l 1 -Boltzmann  distribution. 

3)  The  extremal  focusing  mechanism. 

4)  The  finite  detector  resolution. 

For  the  present  calculations  we  will  Ignore  the  transition  probability  as  a 
function  of  wavelength.  This  variation  is  expected  to  be  minor  with  respect 
to  the  expected  singularities  due  to  extremal  focusing.  The  simulations 
were  done  by  solving  Eq.  (3)  for  a  weighted  distribution  of  several 
thousand  Individual  wavelengths  and  integrating  the  resultant  distribution 
at  each  angle  over  the  finite  detector  aperature  of  3.25  degrees.  These 
calculations,  for  a  beam  of  HD  incident  on  MgO  (001)  along  the  [110]  azimuth, 
are  shown  in  Figs.  47-30.  In  each  figure,  the  upper  section  represents 
scattering  at  a  beam  temperature  of  300  K  and  the  lower  at  80  K.  The 


Intensity  (arbitrary  units) 


Scattering  Angle, 17  (degrees) 


Calculated  extremal  scattering  peaks  for  HD  from  MgO  (001) 
[110]  azimuth.  49.5  degrees  Incident  angle.  (1,0)  series 
of  inelastic  diffraction  peaks.  lop  panel,  300  K  beam; 
lower  panel,  80  K  beam.  Intensity  in  arbitrary  units. 

1\  -  180*-9-9  . 
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Fig.  48.  Calculated  extremal  peaks  for  HC  at  10.67  degrees  Incident 
angle.  [110]  azimuth,  (1,0)  series  of  inelastic  peaks. 
Upper  panel  300  K  beam;  lower  -  80  K.  Intensity  In 
arbitrary  units. 
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50.  (1,0)  extremal  scattering  peaks.  U9. 4*  incident  angle  for 

300  K  and  80  K  beam.  Intensity  in  arbitrary  units. 


99 


calculated  Intensities  are  plotted  vs  "H,  the  experimental  angle.  These 
distributions  were  calculated  as.  sing  an  equal  number  of  molecules  in 
each  initial  rotational  state.  Only  the  initial  velocities  have  been 
Boltzmann  weighted;  the  distribution  of  initial  rotational  states  will 
be  accounted  for  later. 

Consider  for  example,  the  (1,0)  set  of  peaks  for  9q  •  49.5,  shown 

in  Fig.  47.  At  300  K,  the  -t-1,  +2,  and  +3  peaks  all  exhibit  extremal 

scattering.  The  peaks  rise  rapidly  on  the  low  angle  side  and  then  fall 

precipitously  to  zero  at  the  extremal  angle.  These  sharp  peaks  are  what 

remains  of  the  singularity  after  integration  over  the  detector  aperture. 

The  behavior  for  9  ■  10.67*  and  29.7*  shown  in  Figs.  48-49  la  similar. 

Here  the  -1  and  -2  transitions  clearly  show  the  extremal  behavior.  The 

transition  (l,0)-3  is  different.  This  peak  remains  smooth.  We  recall, 

however,  that  the  extremal  singularity  is  associated  with  the  scattering 

of  a  narrow  range  of  wavelengths  near  \  .  For  (l,0)-3  at  this  angle,  K 

m  m 

is  far  from  the  most  probable  wavelength  in  the  beam.  There  are  thus 
relatively  few  beam  molecules  in  the  region  of  the  singularity.  After 
Integration  over  the  detector,  no  evidence  of  it  remains.  An  interesting 
intermediate  case  is  evident  in  Fig.  50  for  “  49.5*.  Here  the  (I , 0) - 1 
shows  a  peculiar  double  peaked  structure:  one  caused  by  the  peak  in  f(X), 
the  other  by  the  singularity  in  .  In  all  cases,  however,  the  extremal 

edge  occurs,  as  expected,  at  9  A ,  where  A  is  the  half  width  of  the 

detector  aperature. 

At  a  beam  temperature  of  80  K,  the  situation  is  similar;  there  are 
fewer  peaks  involved, however.  These  are  sho»v.  in  the  lower  portions  of 
Figs.  47-50.  In  going  to  lower  beam  temperatures,  the  elastic  peaks  have 
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shifted  outward.  Most  of  the  higher  order  rotational  peaks  have  been 
severely  attenuated;  they  have  become  energetically  closed  to  a  signi¬ 
ficant  portion  of  the  beam  population.  Furthermore,  from  Table  1  on 
page  13,  we  see  that  there  is  an  appreciable  population  of  molecules  in 
the  initial  state  only  for  the  -t-1,  -1,  and  -2  transitions.  This  will 
further  attenuate  the  other  transitions.  As  a  result,  the  +  1  loss  peaks 
along  with  the  elastic  peaks  now  clearly  dominate  the  distributions.  There 
is  little  peak  overlap  to  obscure  the  rotational  features.  Note  that 
the  positions  of  the  extremal  edges  do  not  change  as  the  beam  temperature 
is  changed. 

It  should  be  noted  again  that  not  all  rotational  Inelastic  transitions 
are  affected  by  this  extremal  scattering  mechanism.  An  example  of  non- 
extremal  transitions  is  shown  in  Fig.  31,  calculated  for  helium  at  300  K, 
incident  at  i»9.3*.  The  n«^E<0  transitions  are  shown.  The  distributions 
are  all  smooth.  No  evidence  of  singularities  exists. 

From  these  calculations  we  can  make  several  observations  to  guide  us 
in  our  interpretation  of  the  experimental  data.  The  extremal  scattering 
peaks  should  dominate  the  d lstr ibut ion.  This  should  be  particularly  true 
at  300  K,  where  the  peaks  overlap  to  a  significant  extent,  smearing  out  all 
but  the  sharpest  features.  These  calculations  were  done  assuming  equal 
transition  probability  into  each  peak,  except  when  a  channel  was  closed  to 
a  particular  wavelength.  The  disproportionate  height  of  the  extremal  peaks 
is  solely  caused  by  the  extremal  focusing  mechanism;  it  is  not  caused  by 
differences  in  transition  probability.  We  should  be  able  to  resolve  these 
peaks  even  if  only  a  fraction  of  the  incident  molecules  undergo  these 
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Fig.  SI.  Non-extremal  HD  transitions.  The  (1,0)  gain  peaks  and 
(1,0)  loss  paaks  at  49. S  degrees  are  shown.  Extremal 
singularities  are  not  expected  since  4E*n<0.  These 
smooth  shapes  are  typical  of  such  non-extremal  rotational 
transitions.  Intensity  In  arbitrary  units. 


transit  ions .  However,  only  those  transitions  for  which  n-^E>0  will  have 
these  singularities.  Of  these,  only  transitions  having  a  significant 
number  of  molecules  near  \  will  be  strongly  affected. 

The  300  K  HD  scattering  data  shown  in  Figs.  37-41  have  been  labeled 
with  peak  positions  as  calculated  according  to  this  procedure,  on  the 
assumption  of  a  transition  probability  Independent  of  wavelength.  The  peaks 
which  are  expected  to  have  an  extremal  singularity  have  been  labeled  with 
a  sub-s,  to  indicate  we  expect  them  to  be  sharp.  Those  peaks  for  which  the 
extremal  singularity  should  not  be  present  have  been  labeled  with  a  sub-b, 
for  broad.  Only  "s"  peaks  are  clearly  visible  in  the  experimental  distri¬ 
butions.  These  peaks  occur  very  close  to  the  calculated  positions.  The 
strength  of  these  transitions  is  somewhat  exaggerated  due  to  the  peculiar 
shape  of  the  extremal  peaks. 

The  data  for  scattering  at  a  beam  temperature  of  80  K  have  been 
repeated  in  the  upper  portions  of  Figs.  52-34  with  the  extremal  scattering 
calculations  shown  in  the  lower  panels.  We  recall  that  only  the  -el,  -1, 
and  -2  transitions  have  an  appreciable  population  in  their  initial  states. 

Only  the  elastic  and  the  ♦  l  rotational  transitions  are  evident  in  the 
distributions. 

5.3  Scattering  of  H.  and  D. 

The  preceding  data  and  calculations  were  for  HD.  Similar  effects  should 
occur  for  the  scattering  of  tL,  and  D, .  It  is  instructive  to  consider  what 
we  expect  for  end  D,  in  light  of  our  findings  for  HD  scattering. 

Rotational  scattering  is  expected  to  be  stronger  for  HD  for  several  reasons. 

HD  is  a  very  asymmetric  rotor,  resulting  in  stronger  coupling  to  the  surface 
potential.  Furthermore,  as  shown  in  Tables  1-3,  the  energies  of  the  allowed 

J 


J] 


50  60  70  90  90  100  110  120  130  140  150  160  170 

Scattering  Angie,  17  (degrees)  •*.»*, 


Fig.  52.  HD  scattering  from  MgO  (001).  [ 1 10 J  azimuth.  SO  K  beam 

Experimental  data  in  upper  panel.  Calculations  in  lower 
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50  60  70  80  90  100  110  120  130  140  150 

Scattering  Angle, rj  (degrees) 


Fig.  S3.  HD  scattering  at  80  K.  Same  as  Fig.  52  except  49.4* 
incident  angle. 
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rotational  transitions  for  HD  are  considerably  smaller  (255  cal/mole  vs. 

1012  cal/mole  and  512  cal/mole  for  the  lowest  transition)  of  and  D, . 

It  should  be  noted  chat  since  HD  is  an  asymmetric  rotor,  all  rotational 
transitions  are  allowed  by  symmetry.  For  Che  symmetric  molecules  H.,  and 
D^.the  selection  rule  -  even  is  dictated  by  symmetry;  in  other  words, 
there  is  no  ortho-para  conversion. 

The  distribution  of  rotational  states  for  at  300  K  Is  particularly 
unfavorable  for  detecting  rotational  transitions.  65%  of  the  beam  is  in 
Che  J  ■  1  stace, which  can  only  go  to  the  J-3  state,  that  is  a  -2  transition. 
This  requires  1678  cal/mole  transferred  from  translation.  However,  the  mean 
translational  energy  in  the  beam  is  only  1192  cal/mole.  Using  Eq.  (33)  we 
can  calculate  that  at  300  K,  all  rotational  inelastic  diffraction  channels 
are  closed  to  at  least  54%  of  the  beam  molecules.  Only  the  £  1  transitions 
are  energetically  favorable,  but  only  about  107.  of  the  beam  is  in  each  of 
these  states  (J  ■  0  and  J-2). 

The  rotational  distribution  of  D,  at  300  K  is  more  favorable  (Table  3). 
Transitions  up  to  ♦  3  are  energetically  allowed;  the  rotational  distribution 
is  peaked  at  J  ■  2  but  the  other  states  are  also  reasonably  populated.  From 
mere  energy  and  population  considerations,  the  I  -  +1  (J  ■  2  -♦  J  •  0)  transition 
appears  to  be  the  most  favorable. 

Data  for  the  scattering  of  D.,  and  H,  along  the  [  1 10 ]  azimuth  is  pre* 
sented  in  Figs.  55  and  56,  along  with  the  predicted  positions  of  the 
rotational  scattering  peaks,  using  Che  same  procedure  and  notation  as  used 
for  HD.  A  note  on  the  data  presentation  for  these  two  figures  is  in  order. 

For  compactness,  several  scans  at  different  3^  are  superimposed.  Each  is 
shifted  up  from  the  previous  one.  The  baseline  appropriate  to  each  curve 
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Fig.  55.  scattering  from  MgO  (001).  [110]  azimuth. 

300  K  beam,  Only  every  5-th  data  point  marked. 

Each  run  shifted  successively  upward.  Baseline 
indicated  by  horizontal  portion  of  angle  indicating 
incident  angle. 


Fig.  56.  scattering,  as  In  Fig.  50. 
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Is  indicated  by  the  horizontal  portion  of  the  arrow  indicating  the  incident 
angle.  The  full  intensity  of  the  specular  peak  is  indicated  where  the 
peak  extends  off  the  page.  Again,  the  expected  extremal  peaks  have  been 
labeled  with  a  sub-'V. 

We  will  consider  the  D.  data  first.  It  is  immediate ly  evident  that 

4 

there  is  much  less  rotational  Inelastic  scattering.  The  scattering  is 
much  more  like  the  helium  scattering  in  Chapter  4  than  the  HD  scattering 
earlier  in  this  chapter.  There  are  a  few  peaks  which  can  be  identified  as 
extremal  scattering.  The  (2,0) + 1  for  3^  -  49.1*  is  most  clearly  marked. 
Additional  evidence  can  be  seen  for  the  (2,0)-l  for  •  29.0*  and  the 
(l,0)-l  for  3^  ■  9.3*,  but  these  are  very  weak.  As  expected  from  the 
considerations  earlier  in  this  section,  the  -el  rotational  transition  is  the 
most  distinct. 

It  is  instructive  to  compare  this  D ,  data  to  the  helium  data  presented 
earlier.  Deuterium  and  helium  have  the  same  mass  and  thus  the  same  wave¬ 
length.  The  specular  intensity  for  D,  is  lower  than  for  He  by  almost  an 
order  of  magnitude.  This  is  due  to  a  large  extent  to  the  Increased  in¬ 
elastic  scattering  involving  rotational  modes  of  the  molecules. 

The  data  in  Fig.  56  show  even  less  rotational  scattering. 

Little  direct  evidence  of  any  of  the  extremal  transitions  can  be  seen.  The 
Intensity  of  the  specular  peak  is  considerably  more  for  end  D^*  This 
may,  however,  not  all  be  due  to  decreased  Inelastic  scattering.  The  wave¬ 
length  of  H.,  is  v2  longer  than  that  of  or  helium  at  the  same  energy. 
There  are  thus  fewer  elastic  diffraction  channels  open  as  well. 


5.4  Summary 

The  sharp  peaks  which  appear  in  the  scattering  distributions  of  hydro- 
genic  molecules  from  MgO  are  in  keeping  with  the  effects  predicted  by  our 
theory  of  extremal  scattering.  These  peaks  are  actually  singularities  in 
the  differential  scattering  intensity.  For  these  rotationally  inelastic 
transitions,  molecules  with  a  significant  range  of  wavelengths  can  be 
focused  into  a  small  angular  increment.  The  resulting  sharp  peaks  stand 
out  against  the  smooth  background  of  elastic  scattering.  In  a  sense,  the 
use  of  a  beam  with  a  broad  energy  distribution  has,  therefore,  proved  to 
be  an  advantage. 

Extremal  scattering  occurs  for  all  the  hydrogenic  molecules,  but  is 
particularly  strong  for  HD.  It  could,  however,  occur  in  the  scattering 
of  other,  heavier  diatomic  molecules.  Several  things  stand  in  the  way, 
however.  Diffraction,  a  necessary  precondition  for  extremal  scattering, 
has  not  yet  been  observed  for  these  heavier  molecules.  Furthermore,  the 
rotational  levels  of  these  molecules  are  so  closely  spaced  that  the  peaks 
would  not  be  separately  resolvable.  The  occurrence  of  extremal  singularities 
has  been  very  important,  in  allowing  us  to  establish  rotational-trans¬ 
lational  Interchange.  Rotational  inelastic  events  of  this  type  occur  with 
high  probability  for  the  hydrogenic  molecules.  Although  direct  data  are 
not  available  for  other  molecules,  we  expect  rotational  inelastic  scattering, 
as  described  here,  to  play  an  important  part  in  energy  transfer  for  other 


diatomic  molecules  at  surfaces 
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CHAPTER  6 

DIFFUSION  OF  CHEMISORBED  GASES 

The  dynamics  of  chemical  reactions  on  surfaces  is  governed  to  a  large 

degree  oy  the  energetics  of  adsorption,  diffusion,  and  desorptions  of  the 

reactant  species.  In  a  typical  reaction,  adsorbate  molecules  may  have  to 

diffuse  to  an  active  site  or  to  each  other  in  order  to  react.  The  rate  of 

diffusion  relative  to  the  rate  of  desorption,  for  example  ,  may  be  of  critical 

importance.  Little,  however,  is  known  about  the  energetics  of  diffusion  of 

53  54 

even  simple  gases  on  low  index  crystal  surfaces.  ’  One  of  the  problems 
thac  has  stood  in  the  way  of  learning  more  about  diffusion  is  the  need  to 
create  a  spatially  localized  initial  deposit,  the  spreading  of  which  can  be 
studied.  During  the  course  of  this  investigation  of  molecular  beam 
scattering,  we  have  incidently  mastered  this  particular  aspect  of  the 
problem,  and  have  developed  a  quite  general  technique  for  obtaining  data 
on  the  surface  diffusion  of  gases  on  macroscopic  single  crystal  surfaces. 

It  combines  a  molecular  beam  source  with  an  Auger  electron  spectroscopy 
detector.  The  technique  can  be  used  for  any  gaseous  species  which  adsorbs 
readily  on  the  substrate  of  interest,  with  the  notable  exception  of  hydrogen. 
The  substrate,  however,  must  be  in  a  form  to  allow  rapid  heating  and 
cooling.  Although  quantitative  data  are  not  yet  available  because  of 
difficulties  with  some  of  the  instrumentation,  the  proposed  technique  and 
the  equipment  will  be  briefly  described  in  this  chapter. 

/  '  * 

6 . 1  Experimental  Apparatus 

We  wish  to  measure  the  rate  at  which  a  localized  deposit  proceedes 
toward  a  spatial  equilibrium  distribution  by  diffusion  across  the  surface. 

We  thus  require  1)  a  method  of  establishing  a  reproducible  lateral  concen¬ 
tration  gradient  of  the  chemisorbed  gas,  and  2)  a  method  of  measuring  both 
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Che  initial  and  final  concentration  profiles.  We  will  use  our  molecular 
beam  to  create  Che  initial  deposit,  and  a  form  of  scanning  Auger  electron 
spectroscopy  to  measure  it.  To  carry  out  these  measurements  the  scattering 
equipment  in  Figs.  4-5  has  been  modified  as  indicated  by  the  schematic  view 
in  Fig.  57. 

We  are  now  interested  in  the  gas  which  sticks  to  the  surface,  and  not 
that  which  scatters.  We  will  use  Che  molecular  beam  to  form  the  concen¬ 
tration  gradient.  This  is  done  by  slowly  rotating  the  sample  in  an  arc  in 
front  of  the  molecular  beam.  If  done  under  the  proper  conditions,  this 
should  produce  a  narrow  line  of  chemisorbed  gas  running  across  the  crystal. 
The  flux  of  molecules  at  Che  surface  from  the  molecular  beam  is  approximately 
10^  times  the  flux  due  to  the  background  gas,  so  there  should  be  considerable 
contrast  between  the  line  deposit  and  the  background.  The  beam  has  been 
stopped  down  to  approximately  0.020"  diameter  for  these  experiments.  We 
thus  expect  a  line  extending  across  the  crystal  surface  slightly  wider  chan 
0.020". 

Auger  electron  spectroscopy  will  be  used  to  measure  the  lateral  concen¬ 
tration  profile.  In  this  technique,  we  measure  the  energy  distribution  of 
secondary  electrons  when  the  sample  is  bombarded  with  81  3  kV  electrons. 

Each  element  has  a  characteristic  set  of  transitions.  Because  of  the  short 
mean  free  path  of  low  energy  electrons  in  a  solid,  the  technique  is  only 
sensitive  to  the  surface  layer.  Concentrations  down  to  0.1*  can  be 
detected  for  most  elements. 

A  commercially  available  cylindrical  mirror  electrostatic  energy 
analyzer  (CMA)  (PHI,  10-234  G)  is  used  to  detect  the  Auger  electrons.  The 
spatial  resolution  is  limited  by  the  lateral  size  of  the  electron  beam. 
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Our  analyzer  delivers  a  beam  0.004"  in  diameter  at  3  kV.  Scanning  Auger 
electron  spectroscopy  typically  Involves  electrostatically  deflecting 
the  beam  back  and  forth  across  che  surface  with  the  analyzer  fixed.  One 
can  thus  obtain  lateral  concentration  profiles.  The  field  of  view  of  most 
electron  energy  analyzers  is  quite  small,  however;  one  is  typically  limited 
to  areas  less  than  1/2  ram  by  1/2  mm.  That  is  too  small  for  this  experiment, 
as  the  initial  deposit  is  larger  than  that.  We  will  therefore  keep  the 
analyzer  and  the  electron  beam  fixed,  and  accurately  translate  the  sample 
across  the  front  of  the  analyzer,  perpendicular  to  the  line  deposit. 

Although  this  method  is  admittedly  more  cumbersome,  it  has  the  advantages 
thac  line  scans  can  be  done  across  a  large  area  and  that  the  transmission 
of  the  analyzer  will  not  vary  with  position. 

The  experiment  requires  rapid  heating  (to  2400  K)  and  cooling  in  ultra 
high  vacuum.  The  sample  holder  must  be  vibration  free;  also,  it  must  accu- 
rately  translate  in  the  y-dlrectlon  (perpendicular  to  the  line  source)  in 
increments  of  0.001"  for  measurements  of  the  diffusion  profile,  finally,  it 
must  rotate  by  more  than  180*  to  allow  positioning  of  the  sample  in  front  of 
the  Auger  analyzer  or  the  molecular  beam,  as  shown  in  Fig.  57.  A  new  sample 
holder,  incorporating  these  features,  has  been  built.  In  the  present  design, 
the  sample,  in  the  form  of  a  0.010"  thick  wafer,  is  electron  beam  welded  to 
two  0.080"  W  pins.  These  are  firmly  held  in  molybdenum  clamps  which  are  an 
Integral  part  of  the  liquid  nitrogen  cooling  tubes,  as  shown  in  Fig.  58.  The 

sample  is  heated  by  passing  current  directly  through  it.  On  the  order  of  100 

A 

amps  are  required  to  heat  the  crystal  to  2000  K;  \C  cools  from  1400  K  to 
room  temperature  in  »  15  seconds.  The  requirements  of  good  thermal  and 
electrical  contact  as  well  as  rigidity  are  at  odds  with  the  thermal 


Fig.  58. 


Schematic  view  of  sample  holder  for  diffusion  studies 


expansion  of  Che  crystal  when  heaced  Co  2^*00  K.  These  problems  have  pre¬ 
vented  quantitative  experiments  co  date.  Most  difficulties  appear  to  be 
solved,  however.  We  have  been  able  to  create  a  line  deposit  of  molecularly 
adsorbed  v-nitrogen  on  a  W(110)  surface.  As  yet  we  have  not  attempted  to 
deposit  chemisorbed  3-nltrogen,  chat  is,  nitrogen  which  is  held  as  atoms. 
However,  a  brier  indication  of  how  the  experiments  should  be  carried  out 
appears  in  order. 

6.2  Proposed  Experiments 

Measuring  the  diffusion  of  gases  on  a  crystal  involves  the  following 
operations : 

1)  cleaning  the  sample, 

2)  using  the  molecular  beam  to  deposit  a  line  source, 

3)  using  Auger  electron  spectroscopy  to  measure  the  initial 
concentration  profile, 

A)  heating  at  a  fixed  temperature  for  a  fixed  period  of  time, 

S)  measuring  the  final  concentration  profile. 

These  steps  must  be  repeated  for  a  number  of  temperatures  and  time  Intervals 
to  characterize  the  diffusion  process. 

We  are  interested  in  diffusion  in  a  plane,  in  a  direction  perpendicular 
to  a  line  source.  The  initial  concentration  profile  C(x,0),  and  the  final 
concentration  profile,  C(x,t),  are  actually  measured.  Diffusion  proceeds 
according  to  Fick's  second  law,^  which  in  one  dimension  Is 

3t  Sx  (C  ix  ’  * 


(38) 
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where  0  is  Che  diffusion  coefficienc.  We  will  assume  chat  D  is  inde¬ 
pendent  of  concentration;  Eq.(38)  then  reduces  to 


-  D  C  (x,t) 


(39) 


for 


For  an  Initial  distribution  in  the  form  of  a  delta  function,  that  is, 


C(x,c  -  0)  -  Cq5(x)  . 


the  solution  to  Eq.  (39)  can  be  written  by  inspection  as 


C (x , c) 


C 

o 

vZhDc 


exp (-x  / 4Dt ) 


(40) 


The  beam  source  cannot  be  expected  to  create  such  an  initial  distribution. 
For  a  more  realistic  distribution,  we  proceed  as  follows.  By  integrating 
Eq.  (uQ)  over  the  initial  distribution,  we  obtain  the  general  solution 


C(x.t)  -  —2 -  J  f(x+»)  exp(-?2/4Dt)  d5  .  (41) 

V  4r!Dt 

where  f(x)  ■  C(x,t  *>  0)  Is  the  experimentally  measured  initial  concentration 
at  x.  Substitution  into  Flck's  second  law,  Eq.  (39),  confirms  that  this  is 
the  correct  solution.  For  simple  distributions,  this  Integral  can  be  done 
analytically.  Otherwise,  it  is  necessary  to  evaluate  it  numerically.  In 
either  case,  the  experimental  data,  that  is,  C(x,0)  and  C(x,t)  at  a  fixed 
temperature,  are  fit  to  this  solution  using  the  diffusion  constant,  D, 
as  a  variable  parameter. 

Diffusion  on  a  solid  surface  generally  occurs  over  an  activation 
barrier.  The  atom  must  gain  sufficient  energy  to  surmount  the  barrier 
between  binding  sites.  This  is  shown  schematically  in  Fig.  59.  The 
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diffusion  coefficient  can  be  written  in  terms  of  this  activation  barrier, 
* 


as 


-E  */kT 


D(T)  -  DQe  m  ,  (42) 

where  T  is  the  temperature.  From  measurement  of  D(T)  vs  T  we  obtain  the 

activation  barrier  to  diffusion.  Our  diffusing  atoms  thus  provide  us  with 

a  microscopic  probe  of  the  surface  potential. 

The  above  technique  can  be  applied  to  almost  any  gas-solid  system  in 

which  the  gas  chemisorbs  readily  from  the  gas  phase.  The  major  exception 

is  hydrogen,  which  has  no  Auger  transitions.  Some  other  method  would  be 

necessary  to  detect  hydrogen.  We  have  begun  studies  on  the  diffusion  of 

chemisorbed  nitrogen  on  tungsten  (100).  The  adsorption  of  nitrogen  on 

tungsten  has  previously  been  extensively  characterized.  This  makes  it  an 

ideal  system  for  our  initial  studies. 

There  are  numerous  interesting  problems  to  which  this  technique  can 

be  applied.  For  instance,  the  diffusion  of  a  gas  on  various  planes  of 

the  same  material  can  be  studied.  Of  particular  Interest  is  diffusion  on 

highly  anisotropic  planes  such  as  the  W(112),  which  is  composed  of  one- 

56 

dimensional  channels.  Studies  in  the  field  ion  microscope  have  shown 
that  for  metal  adatoms,  diffusion  occurs  only  along  the  channels,  not  over 
the  ridges.  The  correlation  of  diffusion  behavior  of  chemisorbed  gas  with 
substrate  structure  would  provide  us  with  valuable  information  on  microscopic 
surface  forces. 


II 
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